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Summary 


The goals of this research program were to 

(i) determine how microstructural factors, especially the architecture of reinforcing fibers, 
control stiffness, strength, and fatigue life in 3D woven composites; 

(ii) identify mechanisms of failure; 

(iii) model composite stiffness; 

(iv) model notched and unnotched strength; and 

(v) model fatigue life. 

We have examined a total of eleven different angle and orthogonal interlock woven 
composites. Extensive testing has revealed that these 3D woven composites possess an 
extraordinary combination of strength, damage tolerance, and notch insensitivity in compression 
and tension and in monotonic and cyclic loading. In many important regards, 3D woven 
composites far outstrip conventional 2D laminates or stitched laminates. Detailed microscopic 
analysis of damage has led to a comprehensive picture of the essential mechanisms of failure and 
how they are related to the reinforcement geometry. The critical characteristics of the weave 
architecture that promote favorable properties have been identified. Key parameters are tow size 
and the distributions in space and strength of geometrical flaws. The geometrical flaws should be 
regarded as controllable characteristics of the weave in design and manufacture. 

In addressing our goals, the simplest possible models of properties were always sought, in 
a blend of old and new modeling concepts. Nevertheless, certain properties, especially regarding 
damage tolerance, ultimate failure, and the detailed effects of weave architecture, require 
computationally intensive stochastic modeling. We have developed a new model, the “Binary 
Model”, to carry out such tasks in the most efficient manner and with faithful representation of 
crucial mechanisms. 

This is the final report for contract NAS 1-18840. It covers all work from April 1989 up to 
the conclusion of the program in January 1993. 
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Introduction 


Recent advances in resin transfer molding have made possible the manufacture of high 
volume fraction composites from a great variety of dry fiber preforms. The most interesting of 
these have three-dimensional (3D) reinforcement, in which the fibers alone carry loads in each of 
three linearly independent directions. Dry 3D preforms can be manufactured by weaving, braiding, 
stitching, or knitting [e.g., 1-5]. 

3D polymer composites have been pursued in the aerospace industry largely in reaction to 
the vulnerability of conventional laminates to delamination under impact and subsequent failure by 
buckling under in-plane compression. In this regard, they have been successful [3, 4,6,7]. It is 
now becoming clear that they have outstripped this goal by also exhibiting extraordinary properties 
under in-plane loading in their pristine condition. Following their 3D carbon-carbon composite 
ancestors [e.g., 8], certain 3D polymer composites offer large strains to failure 1 , suggesting notch 
insensitivity and damage tolerance far beyond those of 2D laminates. 

In early work under this contract [9], high strains to failure of 3D woven interlock 
composites under uniaxial compression were linked to broad distributions of geometrical flaws in 
strength and space. The geometrical flaws comprise both misalignment of nominally straight load 
bearing tows and certain topological features of the weave architecture. Geometrical flaws lower 
the critical value of the local stress for kink band formation, which was determined to be the 
primary mechanism for failure. 

In tension and bending, damage tolerance is again far in excess of that typical of 
conventional 2D laminates. Consideration of the mechanics of the observed failure events leads to 
general remarks about the origins of such behavior. Once again, geometrical flaws appear to have a 
vital role. 

The mechanisms of damage observed in fatigue are closely related to those seen under 
monotonic loading. Ultimate failure again tends to be associated with broadly distributed damage 
and high energy absorption. 


1 Throughout this report, strain to failure will refer to the strain at which any specimen ceased to bear load when 
tested under stroke control. The strain failure so defined is always equal to or greater than the strain to peak load; it is 
often much greater. 
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All our experimental observations form the basis for models for design and reliability in 
monotonic and cyclic loading. The first generation of these models is also reported here. 

2. Materials 

The subject materials were various 3D interlock weaves impregnated and consolidated with 
epoxy resin. Eleven different composites were studied in all. The weaver’s specifications for the 
reinforcement in each case are listed in Table 1. As also indicated in Table 1, the eleven materials 
can be conveniently divided into two groups according to the degree of fiber compaction achieved 
during processing. 

Table 1 

Weaver’s Specifications for Reinforcement 




Warp Tow Denier 1 


Number of 



Preform 


Weaver 



layers of 

Ends per 

Picks per 

Label 

Architecture 

Material Stuffcrs 

Fillers 

Weavers 

Stuffers 

cm (incli) h 

cm (inch) 0 

(a) Lightly Compacted (Resin: Dow Tactix 138 with 1141 hardenert ) 






/-L-l 

Layer-to- Layer 

AS4 13.8k 

13.8k 

5.9k 

4 

5.1 (13) 

4.4 (11.3) 

i-L-2 

Angle Interlock 

S-2 glass 13.8k 

13.8k 

5.95k 

4 

5.1 (13) 

5.9 (15) 

/- T-l 

Through -the-Thickness 

AS4 13.8k 

13.8k 

5.9k 

4 

4.7 (12) 

5.0 (12.7) 

J-T-2 

Angle Interlock 

S-2 glass 1 3.8k 

13.8k 

5.95k 

4 

5.1 (13) 

5.0 (12.7) 

/- 0 

Orthogonal Interlock 

AS4 13.8k 

13.8k 

5.9k 

4 

4.7 (12) 

5.1 (13) 

(b) Heavily Compacted (Resin: Shell 

RSL-IX95 with HI’ON CURING AGiiNT* W* ) 




A-L-l 

Layer- to- Layer 

AS4 15.8k 

7.9k 

(3.94k,0.66k) d 

4 

5.5(14) 

5.1 (13) 

A- L-2 

Angle Interlock 

AS4 7.9k 

3.94k 

(l.97k,0.6fik) d 

6 

7.1 (18) 

7.9 (20) 

A- T-l 

Throu gh - the-Th ic k ness 

AS4 15. 8k 

7.9k 

3.94k 

4 

5.5(14) 

5.1 (13) 

A-T-2 

Angle Interlock 

AS4 7.9k 

3.94k 

1.97k 

6 

7.1 (18) 

7.9 (20) 

AO-1 

Orthogonal Interlock 

AS4 15.8k 

7.9k 

3.94k 

4 

5.5(14) 

5.1 (13) 

A-O-2 

Orthogonal Interlock 

AS4 7.9k 

3.94k 

1.97k 

6 

7.1 (18) 

7.9 (20) 

a For AS4 fibers, a tow denier of 13.8k corresponds to a filament count (fibers per low) of 21,000, 15.8k lo a count of 24,000, 7.9k to 12,000, 


,0.66k to 1000. 


b number of columns of staffers per cm (inch) in weft direction. 
c number of columns of fillers per cm (inch) in warp direction, 
d The first figure refers to warp weavers, the second to surface warp weavers. 


t Dow Chemical Co., Freeport, Texas 
* Shell Oil Co., Anaheim, California 


All fiber preforms were supplied by Textile Technologies, Inc.t The “lightly compacted 
composites” were processed by the authors and their colleagues by methods described fully in 
Appendix A, using Dow Tactix 138 resin with H41 hardener. The dry preforms in this group 
were relatively loosely woven and minimal compaction pressure was applied to them in the 
through-thickness direction during consolidation. The resulting composites possess relatively low 


t Textile Technologies, Inc., Hatboro, Pennsylvania. 
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total fiber volume fractions. The “heavily compacted composites” were processed by Boeing* 
using Shell 1895 resin with a cure cycle of 30 min. at 149°C and 90 min. at 177°C, followed by a 
further 2 hr. post cure at 177°C [10]. The dry preforms in this group were relatively tightly woven. 
Substantial through-thickness pressure (1.5 MPa) during consolidation produced some further 
compaction and ensured the attainment of predetermined thicknesses for each composite. The 
resulting volume fractions were relatively high. Because of the higher volume fractions and the 
lower weight per unit area of the dry fiber preforms, the heavily compacted composites were only 
half as thick as the lightly compacted ones (Table 2). 

Table 2 


Fiber Volume Fractions 


Composite 

Specimen 

Thickness 

Nominal Volume Fractions 

Total Fiber 
Volume Fraction® 

Fraction by Volume of All Fibers Lying in: 1 
Stuffers Fillers Warp Weavers 


(cm)* 

v,*> 

Vf 

Vw d 

V 

fs 

ff 

fw 

/-L-l 

1,26 

0.14 

0.15 

0.07 

0.35 ± 0.03 

0.385 

0.418 

0.197 

/- L-2 

1.24 

0.14 

0.20 

0.05 

0.370 ±0.0051 0.066 ± 0.004** 

0.347 

0.502 

0.151 

/-T-l 

1.02 

0.16 

0.21 

0.05 

0.466 ± 0.003 

0.381 

0.504 

0.115 

f-T-2 

0.97 

0.18 

0.22 

0.04 

0.408 ± 0.0208 0.044 ± 0.004*> 

0.406 

0.497 

0.097 

1-0 

0.88 

0.18 

0.25 

0.04 

0.483 ± 0.010 

0.387 

0.524 

0.090 

A- L-l 

0.561 

0.38 

0.22 

0.05 

0.620 ± 0.008 

0.587 

0.340 

0.073 

A- L-2 

0.625 

0.33 

0.21 

0.025 

0.557 ± 0.015 

0.580 

0.375 

0.045 

A-T-l 

0.573 

0.37 

0.22 

0.065 

0.61310.003 

0.571 

0.331 

0.098 

A-T-2 

0.577 

0.36 

0.23 

0.035 

0.592 1 0.014 

0.571 

0.369 

0.059 

A-O-l 

0.579 

0.37 

0.22 

0.045 

0.619 1 0.008 

0.586 

0,340 

0.073 

A-O-2 

0.587 

0.35 

0.23 

0.065 

0.593 ± 0.014 

0.545 

0.353 

0.102 


a in direction normal to warp and weft directions. 
b V s 2 volume fraction of stuffer (straight warp) tows. 
c Vf b volume fraction of filler (weft) tows, 
d V w * volume fraction of warp weaver (3D warp) tows. 
e measured by acid digestion. 
f determined from weaver's specifications. 

& graphite fibers. 
h glass fibers. 

The angle and orthogonal interlock architectures are shown schematically in Fig. 1. The 
stuffers (straight warp) and fillers (straight weft) form a 0°/90° array resembling a coarse laminate, 
but with individual tows remaining distinct within each lamina. Figure 1 shows projections on 
planes normal to the filler direction of architectures representative of those in Table 1. Progressing 
down the filler direction, planes containing a full complement of stuffers generally alternate with 
planes containing two warp weavers in the case of angle interlock weaves and one warp weaver in 
the case of orthogonal interlock weaves. The phases of the warp weavers encountered in 
successive planes are varied to minimize short range order. In the thicker orthogonal interlock 


* Boeing Aircraft Co., Seattle, Washington. 
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composites, planes containing a warp weaver are alternated with two layers of stuffers to limit the 
volume fraction of the warp weavers to the desired range. 



q filler (weft) 


stuffer (straight warp) 



warp weaver 


Fig. 1 Schematics of (a) layer-to-layer angle interlock (b) through-the-thickness angle interlock, 
and (c) orthogonal interlock weaves. Numbers indicate typical sequences in which warp 
weavers are encountered on progressing down the filler direction. 
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Table 2 summarizes fiber volume fraction data for all the composites. The nominal volume 
fractions for each kind of tow are deduced from the weaver’s specifications of Table 1 and the 
measured composite thickness of Table 2. Following the ideal geometry of Fig. 1, the stuffers and 
fillers are assumed straight in these estimates, while the warp weavers are assigned a “take up 
factor”, which is the ratio of the measured length of representative warp weavers drawn from a 
woven preform to the length of preform from which they have been drawn. The take up factor is 
typically ~ 1.2 for layer-to-layer angle interlock weaves, ~ 1.3 for through-the-thickness angle 
interlock weaves, and between 3 and 5 for orthogonal interlock weaves. Since the warp weavers 
contribute relatively lightly to the overall fiber content, the precise value of the take up factor has a 
small influence on in-plane property estimates. In practice, some departure from these nominal 
volume fractions is to be expected, because variations may occur in the tightness of the fabric and 
sizing on fibers influences the specified tow denier. Therefore, the total volume fraction, V, 
constituted by all fibers was measured for each composite by weighing the fibers in a known 
volume after acid digestion of the resin (following the procedure in ASTM Standard D3171). The 
value of V differs from the sum of the nominal volume fractions by as much as 10%. Regarding V 
as the more accurate and reliable figure, the volume fractions of fibers in the different kinds of 
tows can be deduced from it with reasonable confidence by assuming that the proportions of all 
fibers in each kind of tow remain as deduced from the weaver’s specifications; i.e., in the ratios 
Vs.VfiVw. The variables f s , ff, and f w in Table 2 denote the fractions by volume of all fibers lying 
in stuffers, fillers, and weavers, respectively, so determined. 

Both the lightly and heavily compacted composites deviate substantially from the nominal 
or ideal pattern of straight stuffers and fillers and approximately sawtooth, square wave, or 
sinusoidal warp weavers. One of the main aberrations in both groups is the distortion of the fillers 
and stuffers near sites where warp weavers wrap around fillers, as in the schematic of Fig. 2. In 
occurrences of this configuration at the specimen surface, the warp weaver is sometimes separated 
from the surface by a pocket of resin, implying that the distortion must have arisen during weaving 
as a result of tension in the warp weaver. In other occurrences, the warp weaver is pressed against 
the surface, implying that the distortion may be the result of through-thickness compaction during 
consolidation. In the lightly compacted composites, instances of the distortion of Fig. 2 tend to 
appear randomly throughout the structure, while in the heavily compacted composites, they tend to 
occur systematically, often at every warp weaver extremum. The magnitude of the distortion varies 
significantly from one weave type to another. 
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SC-0712-C 



Fig. 2 Schematic of the deformation of a stuffer and filler near an extremum of a warp weaver. 

In the lightly compacted composites, other random distortions of stuffers and fillers are 
also found. In the heavily compacted materials, the stuffers are relatively straight away from warp 
weaver extrema. To some extent, this difference is attributable to the difference in compaction 
pressure and to some extent to the fact that the fiber preforms used for the lightly compacted 
specimens were significantly more irregular as received from the weaver than those used in the 
heavily compacted composites. On the other hand, the fillers in heavily compacted materials are 
generally very heavily distorted (Fig. 3). This is perhaps attributable to the fillers being much 
lighter than the stuffers in the heavily compacted materials and to the fact that fillers are more 
directly affected by the distortion of Fig. 2. 

Another significant distortion is crimping of warp weavers, as illustrated in the micrograph 
of Fig. 4. This distortion is generally much more severe in the heavily compacted composites. It is 
obviously the result of through-thickness compaction during consolidation. 

Aberrations from ideal geometry are very important in determining composite properties, 
especially strength, strain to failure, and fatigue life. 
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fillers 


0.5 cm 

Fig. 3 Typical distortions of (a) stuffers and (b) fillers in the heavily compacted composite 
h- L-l. 

(a) heavily compacted T 


tey ' Lilli I I m 



warp weaver 


(b) lightly compacted 



warp weaver 


Fig. 4 Distortion of warp weavers in angle interlock woven composites that are (a) heavily 
compacted and (b) lightly compacted. 
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3. Monotonic Loading 

3.1 Experiments and Observations 

Because our ultimate goal has been to formulate mechanism-based design and reliability 
models, the experiments reported here were not planned to provide comprehensive engineering 
data for all the composites studied. Emphasis lay instead on associating typical macroscopic 
properties with failure mechanisms and identifying microstructural aspects of the composites upon 
which the failure mechanisms depend. Tabulated data usually refer to a single test for each 
composite type in Table 1, the exceptions being those results for which experimental scatter is 
indicated. 

This section reports tests conducted under monotonic loading. Fatigue tests are reported in 
Section 4. Tests comprised uniaxial tension and compression using dogbone specimens, 
compression of short cuboidal specimens between flat platens, and four-point and three-point 
bending of long rectangular bars. The thickness of all specimens was that of the original composite 
panel, with all machining cuts made normal to the original surfaces. Satisfactorily smooth and 
damage-free cut surfaces were obtained using a water jet. 

Dogbone specimens for the lightly compacted composites had straight sided gauge sections 
approximately 2.4 cm long and 1 . 1 cm wide (Fig. 5a). For the heavily compacted composites, 
some experimentation with specimen shape was required to induce failures in the gauge section. 
The final specimen shape, which was almost universally successful in this regard, was that of 
Fig. 5b. Short cuboidal specimens were approximately 2. 1 cm long and 1 cm wide. Bending bars 
were approximately 12.5 cm long and 1 cm wide, with moment arms of 5.7 cm and 2.9 cm for 
three-point and four-point bending respectively. 

In all compression and tension tests, the long axis of the specimen and therefore the load 
axis were aligned with the stuffers.t In bending tests, bending was applied about an axis parallel to 
the fillers. 


t Tests with loading along the filler direction will be reported under our new contract NAS 1-19243. All aspects of 
modeling for this load orientation are qualitatively identical to the modeling reported here. 
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SC-2S37-CS 

(a) (b) 



Fig. 5 Dogbone specimen dimensions for tests of (a) lightly and (b) heavily compacted 

composites. 

3.1.1 Compression Tests 

All tests were conducted under displacement control. For compression tests, the controlled 
displacement was the stroke of the test apparatus. Stress-strain curves are presented in Fig. 6. 

Test data for the first four lightly compacted composites of Table 1 are summarized in 
Tables 3 and 4. These particular materials have high strains to failure under compression, in 
contrast to conventional 2D laminates or stitched laminates, which fail in a brittle manner. The 
ductility is due to the action of geometrical flaws, which consist of either the misalignment of 
nominally straight tows or certain topological features of the reinforcement architecture. Details of 
this crucial idea appear in [ 1 ] and below. 

Test data for compression tests of heavily compacted composites are summarized in 
Table 5. For both dogbone and cuboidal specimens, substantial nonlinearity sets in at less than 
half the peak load or at strains - 0.5%. The figure quoted in Table 5 for the modulus is that for the 
approximately linear regime at low stresses. In the dogbone tests, softening sometimes reflected 
Euler buckling of the entire specimen. While buckling was rarely observed in tests of lightly 
compacted composites, it is much harder to suppress in the heavily compacted composites because 
they are only half as thick. For the cuboidal specimens, initial nonlinearity can be caused by the 
constraint of the specimen ends and Poisson’s effect, which combine to cause the lateral surfaces 
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Stress (MPa) Stress (MPa) 


(a) heavily compacted composites 
cuboidal specimens 



Strain 


(c) lightly compacted composites 
cuboidal specimens 



Strain 


(b) lightly compacted composites 

Hnn-hnne specimens 



Strain 


(d) lightly compacted composites 
dog-bone specimens 



Strain 


Fig. 6 Stress-strain curves for monotonic uniaxial compression in the stuffer direction. 
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Table 3 


Data for Cuboidal 

Test Specimen — Lightly Compacted Composites 


Specimen 

label 


Maximum 

compressive 





V s a 

stress 

(MPa) 

Modulus 

(GPa) 

e m b 

e i/2 c 


//-L-l 

0.28 

230 

35 

0.0095 

0.06 


/-L-2(I) 

0.23 

215 

23 

0.012 

0.09 

Woven 

/- L-2(H) 

0.20 

205 

20 

0.012 

0.12 

composites i 

/-T-l(I) 

0.33 

240 

38 

0.012 

0.16 


/-T- 1(11) 

0.27 

205 

28 

0.012 

0.15 


/-T-i(in) 

0.25 

200 

31 

0.011 

0.16 


1/-L-2 

0.27 

240 

27 

0.010 

0.1 

Tactix 138 / 

e-1 


105 

2.9 

0.079 

d 

epoxy 1 

e-2 


115 

3.06 

0.080 



a Volume fraction of aligned fibers, including straight surface warp weavers. 
•’Compressive strain at maximum load. 

c Compressive strain when post-peak load has fallen to half the maximum . 
d Both tests of epoxy finished with load still exceeding half maximum. 

Table 4 


Data for Dog-bone Test Specimens — Lightly Compacted Composites 




Maximum 




Method of 

Specimen 

Controlled 

compression 

Modulus 

£m a 


strain 

label 

parameter 

stress (MPa) 

(GPa) 

MTl 


measurement 0 

/-L-l(I) 

strain 

210 

41.5 

0.005 


eg 

/-L-l (II) 

load 

210 

28 



si 

/-L-2(I) 

load 

140 

25 

0.006 


c g 

/-L-2(II) 

strain 

155 

35 

0.005 

0.05 

eg 

/-L-2(in) 

strain 

160 





/- T-l(I) 

strain 

165 

32.5 

0.005 

0.14 

eg 

/- T-l(II) 

load 

150 




eg 

/-T-2(I) 

strain 

205 

35.5 

0.006 

0.03 

c g 

/-T-2(II) 

strain 

195 

38.5 

0.005 

0.04 

si 


Compressive strain at maximum load. 

Compressive strain when post-peak load has fallen to half of maximum. 

C “cg” indicates a clip gauge, “si” steroimaging; in both cases, using a 1.2 cm gauge length. 
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Table 5 


Compression Tests of Heavily Compacted Specimens 


Maximum 

Compressive Modulus 

Preform Label Stress (GPa) e e* 

(MPa) 


Dogbone Specimens 


fi-L-1 

670 

88 

0.008 

— 

h-L-2 

695 

81 

0.009 

— 

h- T-l 

570 

83 

0.0078 

— 

h-T-2 

540 

72 

0.0102 

— 

h- 0-1 

635 

72 

0.0107 

— 

h- 0-2*> 

555 ± 15 

70 ±5 

0.008 

— 

Cuboidal 

Specimens 




f i-L-1 

470 

61 

0.0108 

0.0143 

h- L-2 C 

— 

— 

— 

-- 

/ i-T-1 

515 

52 

0.0112 

0.0146 

h-T-2 

440 

42 

0.0125 

0.025 

h- 0-1 

505 

52 

0.0142 

0.0285 

h- 0-2 c 

— 

— 

— 

-- 


a Strain at which the load falls to half the peak load, 
b Two tests, both of which were terminated by Euler buckling. 
c Not tested. 


to barrel outwards, reducing stiffness because of the concomitant curvature of stuffers. 
Nonlinearity at strains much above 0.5% is also caused by microcracking. 

The strengths of the heavily compacted specimens in compression are two to three times 
those reported for lightly compacted specimens (Tables 3-5). Part of the increase is attributable to 
the higher fiber volume fractions achieved in heavy compaction; but other factors are also 
significant, as discussed below. 

On the other hand, the strains to failure of the heavily compacted specimens are 
substantially lower than those of the lightly compacted materials. The dogbone tests of heavily 
compacted materials reveal essentially brittle behavior, with negligible loads recorded at strains 
above that of the peak load. However, the apparent brittleness may have been exaggerated by 
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failure of the test apparatus to maintain stable stroke control after specimen buckling. The cuboidal 
specimens carried significant loads at strains up to ~ 3%, which far exceed the failure strains of 
conventional or stitched 2D laminates, but remain much lower than those recorded for the lightly 
compacted composites. 

In contrast to lightly compacted materials, the heavily compacted materials are prone to 
delamination failure. While partial delaminations were previously observed in dogbone tests of the 
lightly compacted materials, the associated opening displacements were smaller and they were not 
the mechanism of ultimate failure. Ultimate failure of the lightly compacted materials was always 
the result of kink band formation in individual stuffers and usually occurred at much higher strains 
than delaminations. Figure 7a shows typical damage at high strain (~ 15%) following widespread 
kinking in a cuboidal specimen of the lightly compacted composite Z-T-l. In dog-bone tests of the 
heavily compacted composites, ultimate failure often occurred so quickly after delamination that it 
was difficult to separate the two events. In such cases, the delamination itself remained one of the 
primary manifestations of damage. Figure 7b shows a failed cuboidal specimen of the heavily 
compacted composite h-L-l, where damage consists essentially of brooming of stuffers at one end 
of the specimen, with concomitant lateral displacements accommodated down the specimen by 
delaminations. In this particular example, it is not evident that kink bands ever formed. In most 
other tests where extensive delamination occurred, however, kink bands appeared to be triggered 
by delaminations and contributed significantly to failure. 

Delamination and other factors influencing strain to failure will be discussed further in 
Section 3.2. 
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(a) 


2 mm 
(b) 




Fig. 7 Comparison of compressive failure of cuboidal specimens of (a) a lightly compacted 
composite (/-T-l) and (b) a heavily compacted composite (/i-L-1) that failed by 
delaminating. 

3.1.2 Tension Tests 

Uniaxial tension test data are summarized in Table 6 and representative load-strain histories 
are shown in Fig. 8. The strains of Table 6 and Fig. 8 were measured by a clip gauge of 1 cm 
gauge length attached to one of the machined sides of the specimen. The output of the clip gauge 
was the controlled displacement during the test. Nonlinearity is evident at stresses over half the 
peak load for both lightly and heavily compacted composites. However, whereas nonlinearity 
sometimes appears to be reversible, elastic deformation in compression, in tensile loading it is 
associated with resin microcracking. The heavily compacted materials are far stronger than the 
lightly compacted materials, with strengths often exceeding 1 GPa. The strains to failure are 
ostensibly large, much higher than those of 2D laminates, but the quoted numbers and the 
abscissae of Fig. 8 should be interpreted with care. They do not reflect material constants, but are 
strongly influenced by whether large tensile cracks fall within the clip gauge field, as discussed 
below. 
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Table 6 


Tension Test Data — Loading in Staffer Direction 


Preform Label 

Peak 

Load 

(MPa) 

Measured 

Modulus 

(GPa) 

Strain to 
Failure 3 

Rule of Mixtures 
Modulus b 
E x (GPa) 

/-L-l 

300 ±50 

30 ±6 

0.015 ± 0.005 

36 

/- L-2 

240 ± 10 

28.5 

0.03 ± 0.01 

41 

/-T-l 

350 

27 

0.04 

48 

/- T-2 

240 

39 

0.02 

49 

1-0 

390 ± 30 

30 ±2 

0.022 ± 0.002 

51 

h- L-l 

980 ± 20 

85 ±8 

0.018 

96 

} i-L-2 

935 

80 

0.039 

85 

h- T-l 

840 

79 

0.038 

92 

h- T-2 

895 

72 

0.013 

89 

/i-O-l 

1070 

88 

0.013 

96 

h- 0-2 

850 

69 ± 5 C 

0.013 ± 0.001 

86 


a Interpret with caution! See text. 
b From Eq. (3), with A s = f s . 

c Bilinear stress-strain curves: the figure shown for the modulus is the average over the whole 
range of the test. 


Detailed examinations of tensile failure mechanisms were made for all composite types. 
Photographs of gauge sections taken during tests were examined by stereoscopy to reveal relative 
displacements across cracks, strain fields in intact regions, and overall pictures of damage. The 
principal damage events in both the lightly and heavily compacted composites are matrix cracking 
(both tensile and delamination), tow rupture, and tow pull-out. 

The main load drops observed in stress-strain records, e.g.. Fig. 8, were caused by the 
rupture of one or more aligned tows. When a tow ruptures, it almost always fails across its entire 
cross-section. Matrix cracking around the circumference of the failed tow then debonds it from the 
surrounding composite, so that any stress concentration is minimized and neighboring aligned 
tows commonly remain intact. Sliding along circumferential debond cracks typically extends 
several mm from the location of the rupture. 

The rupture of an aligned tow is also associated with matrix cracking that separates fillers, 
which are the orthogonally disposed tows in the current experiments (Fig. 9). Inspection of 
specimen sections reveals that such interfiller cracks begin at strains exceeding ~ 0.6% for lightly 
compacted composites and ~ 1% for heavily compacted composites. They become widespread after 
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SC-305M 




Fig. 8 Stress-strain records for dogbone specimens of (a) lightly compacted and (b) heavily 
compacted composites under unaxial tension. 


loading to high strains in all specimens studied, with the unexplained exception of the lightly 
compacted material l-O. The layers of resin between fillers are clearly much weaker than the fillers 
themselves, since the fillers are rarely seen to fail internally. The interfiller cracks are analogous to 
the multiple cracks found in the 90° plies of 0790° laminates, except that their spacing is dictated by 
the filler size rather than the mechanics of stress relief. 
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Fig. 9 The separation of fillers under tension in the stuffer direction in an orthogonal interlock 
composite (ft-O-1). Arrow at A shows simple fracture between fillers. Arrows at B show 
instances of fracture followed by falling away of pieces of resin. 

As the number of ruptured stuffers increases, interfiller cracks on one or more planes 
develop large openings, until a macroscopic “tension crack” is evident (Fig. 10). Such tension 
cracks may traverse the whole specimen, but since stuffer failures are not generally coplanar, they 
remain bridged by intact stuffers. One tension crack will commonly cause ultimate failure when the 
bridging stuffers fail and are pulled out of the fracture surfaces (Fig. 10). Tension cracks that have 
been caught within the field of view of photographic records may have opening displacements of 
several tenths of a mm when the load is still near the peak load. Pullout distances for bridging 
stuffers may exceed 1 cm. The case of Fig. 10b shows broken tows that sustained weak pullout 
loads when the total displacement equalled the length of the specimen gauge section. 

Tension cracks do not always cross the whole specimen. When viewed on a cut side of the 
specimen, they are seen occasionally to terminate at a delamination crack running parallel to the 
load axis between a layer of stuffers and a layer of fillers. Ultimate failure may then consist of 
separation of the specimen along a path comprising the first tension crack, the delamination crack, 
and a second tension crack traversing the rest of the specimen. The two tension cracks may be 
offset from one another by ~ 1 cm. More complex systems of tension and delamination cracks 
might be expected in large specimens and are common in bending (see below). Reciprocally, when 
the failure path involves delamination and offset tension cracks, specimens under uniaxial loads 
develop local bending because of broken symmetry. 

The strains reported in Table 6 and Fig. 8 will obviously depend strongly on whether or 
not a tension crack falls within the field of the clip gauge extensometer. If one does, then high 



distinct fillers 
separated by tension crack 


scp. 064 «e.a 2 ii»* 


(a) 


bridging 

•tuffer 



opening : : 5 mm 
of tension 


crack 


load 


axis 



Fig. 10 (a) A lightly compacted composite (/- T- 1 ) and (b) a heavily compacted composite (h-O- 1 ) 
after testing in uniaxial tension. 

strains to failure will be recorded. If not, the measured strain may actually decrease when a tension 
crack forms elsewhere and relieves the load. In this case, truncated stress-strain curves result, such 
as those for materials h-L-2 and h-O in Fig. 8(b). Thus the strains to failure quoted in Table 6 
should not be regarded as material constants, but merely suggestive of high damage tolerance. A 
meaningful quantification of damage tolerance should be based on laws for tow pullout (vide 
infra). 


While stuffers account for most of the composite strength in the orientation chosen for the 
experiments, significant bridging of tension cracks may also be supplied by warp weavers. 
Because they are initially curved, warp weavers qua bridging ligaments are relatively soft, but 
survive high crack opening displacements. 
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Matrix cracks form not only normal to the applied load, but also to a lesser extent at other 
angles under the influence of the local orientation of tows. Significant matrix cracking is usually 
found prior to any stuffer failures, typically at loads above half the peak load, and is correlated 
with macroscopic softening. 

There is some evidence of residual stresses in the matrix. Sensitive strain measurements* 
reveal net compressive strains in the matrix in unloaded cracked specimens relative to the same 
specimens when unloaded and undamaged, which suggests that the cracks have allowed partial 
relaxation of tensile residual stresses. The relaxation of residual stresses is also thought to be 
responsible for the observed normal opening of delamination cracks that lie parallel to the load 
axis. 

3.1.3 Bending Tests 

In the presentation of bending test data, all stresses and strains refer to those that would 
exist in the outer layer of material on the tensile side of the specimen if the material was 
homogeneous, given the load and load point displacement. Some simple remarks on the effects of 
inhomogeneity on strength follow in Section 3.2. 

Figure 1 1(a) shows a four-point bending test of a lightly compacted layer-to-layer angle 
interlock specimen. Failure began at the center of the beam with kink band formation on the 
compressive side followed by tow rupture on the tensile side. Very limited delamination occurred 
on the tensile side only. The material exhibits considerable ductility. Figure 11(b) shows three- 
point bending data for the same material. As indicated in the figure, the main load drops are 
correlated with the kinking and rupture of stuffers. Once again delamination is very limited. 

Four-point bending tests of some heavily compacted materials were spoiled by extensive 
delamination between the outer loading pins. Since shear was evidently a principal failure mode for 
these composites, all further tests were conducted in three-point bending. With some variance from 
test to test duly noted, representative load-strain records are shown in Fig. 12, where significant 
load drops are annotated with corresponding observations. Of the heavily compacted materials, 
both the layer-to-layer angle interlock (Fig. 12a) and orthogonal interlock (Fig. 12b) types showed 
significant delamination, although it was less in the latter. Both types of through-the-thickness 
angle interlock (Figs. 12c and 1 2d) showed no significant delamination. Where delamination 
occurred, it frequently ran from the center of the beam all the way to one end and was always 

* Stereoscopy and automated differential displacement analysis. See [II] and [12] for details of these methods. 
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Strain 

Fig. 1 1 Test results for lightly compacted layer-to-layer angle interlock specimens (type /-L-l of 
Table 1) in (a) four-point bending and (b) three-point bending. 

accompanied by a large load drop. Shear displacement discontinuities of up to 1 mm are commonly 
observed across delaminations that have run the specimen, yet the delaminations remain bridged by 
intact warp weavers (Fig. 13). Some striking periodic patterns were produced on the upper and 
lower surfaces of the specimen in such cases, where the near-surface extremities of bridging warp 
weavers have been pulled down strongly into the body of the composite (Fig. 13). 

In summary, bending tests reveal much the same failure mechanisms for stuffers as 
uniaxial tests, namely kink band formation on the compressive side followed by tow rupture on the 
tensile side. In correspondence with behavior under unaxial compression, the strain to failure in 
bending tends to be larger for the lightly compacted composites than for the heavily compacted 
composites. Substantial delamination does not occur in lightly compacted composites in bending. 
Whether it occurs in heavily compacted composites depends on the reinforcement architecture. 
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Stress (MPa) Stress (MPa) 


delamination 


I 1 1 1 1 1 1 1 

kink band 


complex failure 
including limited 


.01 0.015 

Strain 


- kink band 
failure 


0.01 

Strain 


Fig. 12 Test results for various heavily compacted composites in three-point bending. 




Surface exposed by axial 
displacement of stuffer (-1 mm) 




delamination crack 


Fig. 13 A delamination crack created in a specimen of type h-O-l by three point bending, with 
large mode II crack opening displacement. The upper surface shows systematic damage 
caused by warp weavers that bridge the delamination crack. ' ° 













3.1.4 Compression After Impact 

Tests of compressive strength after impact were made for some of the heavily compacted 
composites. Rectangular specimens were impacted by Professor George Springer and colleagues at 
Stanford University. The impact was effected by a Teflon projectile fired by an air gun. Several 
velocities were used and the damage produced probed by x-ray radiography. A circular zone of 
damage of about 1 cm in diameter was found for the highest velocity, which corresponded to an 
impact energy of ~ 10 ft-lb (13.5 Nt-m). 

Table 7 compares residual strengths after impact with the data for pristine dog-bone 
specimens. Ignoring the case that terminated prematurely via Euler buckling, the loss of strength is 
~ 20%. This is close to the reduction that would be expected from net section considerations. The 
impact damage affects stuffers within a volume of width approximately 1 cm and depth 
approximately 2 mm, while the impacted specimens were approximately 2.5 cm wide and 0.5 cm 
thick. If all stuffers within the impact zone are assumed to have zero strength and stiffness, the 
average load on the remaining stuffers must rise by a factor 1.25/1.05 ~ 1-2. Thus the strength 
should fall by - 20%. Further justification for using net section estimates appears in 
Section 3.2.3. 


Table 7 

Compression After Impact 


Composite 

Label 


Pristine Strength 3 
(MPa) 


Compressive Strength 
After Impact b 


h- L-l 670 345 c 

h- T-l 570 510 

h-T-2 540 435 

h- O-l 635 470 


a From Table 5 - dog-bone specimens. 

b 10 ft-lb (13.5 Nt-m) applied by Teflon projectile. 

c Specimen failed by Euler buckling. 
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3.2 Theoretical Estimates of Properties 
3.2.1 Modulus 


Reasonable estimates of Young’s modulus in the loading direction (i.e., the direction of the 
stuffers) can be found by appropriate rules of mixtures. For this purpose, it is expedient to separate 
fibers into those lying in stuffers and all others, not discriminating fillers and warp weavers among 
the latter. Let V be the volume fraction constituted by all fibers and f s the fraction of all fibers that 
lie in stuffers. If A s is the area fraction occupied by stuffers on a section normal to the load axis, 
then the axial Young’s modulus, E s , for a stuffer is by rule of mixtures 

(l) 


where Ef is the axial modulus of the fibers and E r the modulus of the resin, which is assumed 
isotropic. The remaining resin and fibers, i.e., those in the fillers and warp weavers, may be 
considered to form an effective medium whose modulus, E m , may be approximated in the direction 
parallel to the stuffers by* 


E 


m ~ 


(i-fs)v i , 

. 1-A S E f 


i (l-fs)V ) 1 I ' 1 
1-A S IeJ • 


( 2 ) 


Young’s modulus, Ej, of the entire composite in the direction of the stuffers is approximately 


E, ~ A S E S + ( 1 - A s ) E m . (3) 

If Ef » E r , as for graphite/epoxy composites, the dominant term in Ej is thus f s VEf, which is 
independent of A s . This has practical importance, because it is often difficult to assign an area 
fraction to stuffers when examining cross-sections of composites. The cross-sections of stuffers 
and indeed all tows are often heavily distorted and variable; and measuring A s is further 
complicated by resin-rich regions between tows. On the other hand, V can be measured accurately 
by acid digestion, while f s is known accurately from the records of the weaver. Since Ei is weakly 
dependent on A s /f s for feasible values of A s , it is expedient here to use the approximation A s = f s . 


* Equation (2) is a reasonable approximation for fillers, which are loaded transversely by loads parallel to the 
stuffers. Warp weavers, on the other hand, lie roughly at 45° to the stuffers and might be thought to offer stiffer 
resistance than intimated by Eq. (2). However, warp weavers constitute a relatively small volume fraction and their 
effective modulus is often reduced by crimping. Thus the error should be minor. 
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Using the manufacturer’s datum Ef = 250 GPa for AS4 fibers, the values Er — 3.5 GPa 
for the Shell 1895 resin and E r = 2.8 GPa for Taxtix 138 with H41 hardener (both measured by the 
authors), and the fiber volume fractions of Table 2 with A s = f s , one finds from Eq. (3) the 
estimates of E x for loading in the stuffer direction that are shown in the right-hand column of 
Table 6. (In estimating Ej for composites /-L-2 and /-T-2, the significance of having S-2 glass 
rather AS4 warp weavers was assayed by letting the modulus for the S-2 glass warp weavers be 
either zero or equal to that of AS4 fibers. The values so computed for Ei differed by less than 1%; 
which one would expect since the warp weavers have low volume fraction and influence Ei via the 
relatively small term Em only.) The estimates are 20-50% higher than the experimentally measured 
values for the lightly compacted composites and 10-20% higher for the heavily compacted 
composites. These discrepancies are consistent with rough estimates of the effects of the stuffer 
waviness noted in Section 2. Waviness leads to lower effective moduli, because the reaction of a 
wavy tow to an aligned load is not merely to extend but also to straighten. 

3.2.2 Strength 

Compression 

Kink band formation in fibrous polymer composites is mediated by shear deformation of 
the matrix [13,14]. Under the simplest assumptions concerning the geometry of the kink band, the 
critical axial stress for kinking, Ok, is given by [13-15] 

a k = V4> , (4) 

where t r is the critical stress for shear flow in the resin and <|> is the angle by which the fibers are 
misaligned with respect to the applied load axis. Equation (4) obtains for an infinite array of 
uniformly misaligned fibers, whereas the kink bands of interest here exist within single, finite tows 
whose misalignment varies continuously and randomly along their length. The extent to which 
such details influence the criterion for kinking has not been studied. With such doubts set aside in 
the spirit of identifying trends, some useful numbers emerge when measured values of x T and <(> are 

substituted into Eq. (4). 

The critical shear stress, x r , was deduced to be ~ 75 MPa from tension tests on ±45 
laminates (see Appendix B). The misalignment angle <j> was measured on photographs of stuffers 
exposed on several sections cut perpendicular to the filler direction, since the most severe stuffer 
deflections always occur normal to fillers. The sections revealed all the stuffers in a volume of 
material approximately equal to that of the gauge section in dog-bone tests (as in Fig. 5). Since 
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there are rarely more than two load drops prior to peak load in compression tests, the attainment of 
peak load was somewhat arbitrarily equated to the occurrence of the second kink band. Thus the 
angle <|> corresponding to peak load was taken to be the second greatest misalignment angle 
measured anywhere in a given set of sections. The results of these procedures are shown in 

Table 8. The critical value, o*, k) , of the applied stress at which kink bands first form is related to 
a k by 


Y7 

<4 k) =£r-<*k = A s a k , (5) 

where A s is the area fraction constituted by stuffers on planes normal to the stuffer direction, which 
can be approximated by the fiber fraction f s of Table 2. The values of o* k) are encouragingly 
consistent with the experimental results of [9] and Tables 3 and 4, reinforcing the view that kink 
band formation is the essential mechanism of failure. 


Table 8 


Estimates of First Kink Band Formation 


Material 

T r 

(MPa) 

♦ 

(radians) 

a k = V<J> 

(MPa) 

As 

<4 k) = A s a k 
(MPa) 

Measured 

Ultimate 

Strength 3 

(MPa) 

/-L-l b 

75 

. 





/-L-2 

75 

0.21 

380 

0.385 

140 

150±10 

/-T-l 

75 

0.23 

350 

0.41 

135 

155±5 

/-T-2 

75 

0.17 

460 

0.38 

165 

200±5 

1-0 

75 

0.26 

300 

0.39 

110 

_b 

h-L-l 

75 

0.052 

1520 

0.58 

850 

670 

h- L-2 

75 

0.087 

915 

0.58 

495 

695 

h-T-l 

75 

0.087 

915 

0.57 

485 

570 

h- T-2 

75 

0.087 

915 

0.59 

505 

540 

h-O-l 

75 

0.070 

1145 

0.59 

635 

635 

h- 0-2 

75 

0.070 

1145 

0.54 

575 

555±15 c 


a From Tables 4 and 5. 
b Not measured. 

c Experiment gives a lower bound for c4 k) , since failure was by Euler buckling. 
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Tension 


The manufacturer quotes a strength of 4 GPa for AS4 fibers.* The stress, O s , in the fibers 
in the aligned tows (stuffers) is approximately related to the applied stress a a by 



with Ei given by Eq. (3). For the values given above for Ef and the estimates of Table 6 for Ej, 
Eq. (6) would imply that tow rupture should occur at tensile stresses of 500-700 MPa for the 
lightly compacted composites and ~ 1.5 GPa for the heavily compacted composites. In fact, the 
measured composite tensile strengths are roughly half to two-thirds of these estimates (Table 6). 
There are several likely causes of strength degradation. 1) Prior to weaving, AS4 yarns are 
wrapped with polyvinyl alcohol threads to prevent spreading. The polyvinyl alcohol is removed 
with hot, distilled water when the weaving is completed. Some chemical degradation of the fibers 
is possible. 2) The rigors of the weaving process are always likely to damage tows. 3) As shown 
in Section 2, both the weaving and consolidation processes result in deformation of tows. 
Nominally aligned tows that are not straight will have a lower effective modulus, leading to 
overloading of neighboring tows that are straight. 4) A misaligned, kinked, or nonuniformly 
squashed tow may be weaker than a straight tow. 5) Lateral loads are induced on aligned tows by 
the reaction of warp weavers to in-plane tension, which may reduce their strength. 

Some of these factors also differentiate 2D laminates of plain weave fabric and 
conventional, nonwoven 2D laminates. Prior studies of graphite/epoxy composites in these classes 
have shown that the strength of 2D woven laminates under loads aligned with the warp tows is 
15-25% lower than in comparable 0790° laminates, while differences in fiber volume fraction are 
only ~ 10% [16]. Further research is required to determine which factors account for the strength 
reduction in 3D woven composites being considerably greater than this. For the moment it is 
enough to say that existing 3D woven composites are probably far from optimal in tensile strength. 

Bending 

The ratio of the average stress a s , in the outermost layer of stuffers in pure bending to the 
nominal stress, a surf , quoted in Figs. 10 and 1 1 is approximately 


* Data sheets of Hercules, Inc., Wilmington, Delaware. 
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( 7 ) 


_Es_ = i -2^-^ , 

(jsurf t t 


where t s is the thickness of a stuffer ply, tf is the thickness of a filler ply, and t is the specimen 
thickness. The ratio of the ply thicknesses is approximately 


tf _ n ff 
t s n + 1 f s 


( 8 ) 


while 


n t s + (n + 1) t f = t , (9) 

where n is the number of layers of stuffers. For the values of f s and ff in Table 2, Eqs. (7) through 
(9) yield a s /o surf ~ 0.7 for the composites for which n = 4 and 0.8 when n = 6. 

Multiplying the peak loads in Figs. 1 1 and 12 by these values of the ratio Oj/a 8 ^ provides 
estimates of the local compressive stresses at which kink bands formed. The resulting stresses, 
shown in Table 9, are remarkably higher than the ultimate strength under uniaxial compression in 
all cases (Tables 3-5). 


Table 9 


Critical Compressive Loads in Bending Tests 


Critical Value of Nominal Critical Value of Estimated 

Surface Stress Stress in Outermost Stuffer 3 

Preform Label (MPa) (MPa) 


/-L-l 
h-L-l 
h-T-l 
h-T-2 
h- 0-1 


285 ± 25 
1035 ±35 
1070 
895 
930 


200 ±20 
725 ± 25 
750 
625 
650 


3 from Eq. (7) 


Similarly high compressive strains in unidirectional graphite/epoxy composites under 
bending have recently been reported elsewhere [17]. In [17], the possibility was considered that 
the high compressive strains might be enabled by constraint of microbuckling by the curvature 
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induced by bending. This mechanism was ultimately rejected in [17] on the grounds that 
microbuckling should occur in the plane of the surface, rather than out of the surface; it would be 
unaffected by the curvature of the specimen and would therefore occur at the same lower stress that 
would be expected for uniaxial loading. However, the universal observation for the 3D woven 
composites studied here is that microbuckling of near-surface stuffers does involve deflections out 
of the surface. For the orientation in which specimens were tested here in bending, in-plane 
deflections of stuffers on the compressive surface of the specimen would have to occur in the weft 
direction and would therefore be strongly resisted by fillers. 

The simple expression of Eq. (4) for kink band formation was derived for an infinite array 
of fibers. No analogous criterion is available to account for the presence of a free surface. 
However, one might conjecture that buckling of a stuffer out of a surface, a precursor to kinking, 
will occur most readily at locations where tow waviness lends the stuffer maximum convex 
curvature when viewed from outside the specimen. This conjecture is supported by moir<§ 
interferometry of woven composites in uniaxial compression. Such figures show segments of 
stuffers arching out of one specimen’s surface at locations where kinking subsequently occurred. 
In bending, the curvature of any tow segment will become the sum of that associated with tow 
waviness and the concave curvature produced on the compressive surface by the applied load. The 
maximum convex curvature in a given specimen will be reduced. The stress required to force some 
segment of a tow to arch out of the specimen will thus be raised; and if such arching is a 
prerequisite for kinking, the critical stress will also be raised. Further work is required to test 

this conjecture. 

3.2.3 Ductility and Notch Sensitivity 

Since it dealt with the failure of individual aligned tows in the stress field of an undamaged 
composite, the discussion of the previous section pertains strictly to the critical stress, in either 
compression or tension, for isolated tow failures. How the material progresses thence to ultimate 
failure depends on the distribution of flaws and how loads are redistributed around the site of each 
local failure event. The possibility of composite strength exceeding the stress of first tow failure, 
high strain to ultimate failure, and notch insensitivity are all favored by flaws that are widely 
distributed both spatially and in strength and by mechanisms of stress relaxation around each tow 
failure site. 

In compression, the essential flaws are geometrical in nature, comprising random 
misalignments, which are reflected in <|> in Eq. (4), and certain topological characteristics of the 
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reinforcement, which cause lateral loads in aligned tows. Lateral loads knock down the critical 
stress for kink band formation [14,15]. Kink bands are most commonly observed where 
misalignment is greatest and often near configurations such as that depicted in Fig. 2. (Whether the 
correlation of kink bands with the deformity of Fig. 2 is attributable to lateral loads or whether it 
arises wholly because tows are exceptionally misaligned there is impossible to tell from the 
experimental observations made so far.) 

In tension, the nature of the essential flaws has not been established. Geometrical 
irregularity may be a factor, as surmised in the preceding section, but the mechanics of such effects 
have not been studied. Whereas kink band formation in compression is not influenced to first order 
by the inherent properties of the fibers [e.g., 13-15], that is not so of strength in tension. Inherent 
fiber variance will lead by itself to a Weibull distribution of tow strengths. 

Stress redistribution around a site of tow failure has now been observed in both 
compression and tension to be moderated by debonding of the failed tow from the surrounding 
composite. Accurate strain mapping by stereoscopy reveals large discontinuities in the axial 
displacement across the debond. It is likely that the implied sliding is governed by friction. The 
observed sliding lengths are typically many times the tow diameter, which invites analysis by the 
usual, simple shear lag theory. Let Tf denote the critical shear stress at which sliding is activated. 
Assuming Tf to be uniform over the debond and that the kind band supports negligible axial load, 
the sliding length, / s , satisfies 

s/ s Xf=a t A , (jo) 

where s and A are the circumference and cross sectional area of the tow and a t is the axial stress it 

bears beyond the sliding zone. The degree of stress concentration in tows neighboring a failed tow 
diminishes with increasing sliding length and therefore increases with the magnitude of tf . In a 3D 
composite, Tf will be strongly influenced by triaxial local stresses generated even under uniaxial 
loading by the through-thickness reinforcement (warp weavers). For fixed Tf, the sliding length 
also increases with the tow denier, which is reflected in Eq. (6) by the ratio A/s. The fact that large 
slip distances are observed in the subject 3D composites is thus a direct consequence of the coarse 
scale of the reinforcement. 

Tension 

In the tensile tests reported above, ultimate failure was associated with a band of damage 
traversing the specimen. The band can be viewed as a dominant crack bridged by tows being 
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gradually pulled out of its fracture surfaces (e.g.. Figs. 9a and 9b). In all the tests reported here, 
this crack spanned the specimen as soon as it was observable. However, for a sufficiently large 
specimen, one might envisage such a bridged crack emanating a limited distance from a stress 
concentrator such as a hole. The response of the material will then depend on the nature of the 
bridging or cohesive zone. The detailed mechanics of the failure process are complex and 
computationally challenging, but the essence will be as schematized in Fig. 14. 
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Fie. 14 (a) Conjectured distribution of damage in a large specimen under tension. 

(b) Conjectured tractions acting over central damage plane vs crack opening displacement 
(or half the displacement discontinuity across the damage plane). 
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Figure 14a shows a process zone consisting broadly of two parts. Furthest from the stress 
concentration, damage will consist of interfiller cracks between intact stuffers, a “zone of matrix 
cracking”. Nearer the stress concentrator, stuffers will have failed. Since the stuffer failures will 
not all coincide with the main fracture plane, a “zone of tow pullout” will exist, within which failed 
stuffers continue to transfer loads via friction. When the damage has propagated far enough, the 
stuffers nearest the stress concentrator will pull entirely out from one fracture surface, forming a 
traction free zone. 

The fracture mechanics of this crack system are determined by the relation between the 
tractions p acting across the crack and the crack opening displacement u (or displacement 
discontinuity 2u). The problem is analogous to that of the pullout of fibers of random strength 
from a brittle matrix composite [18,19]. Qualitative analysis leads to a relation p(u) schematized by 
the dashed curve of Fig. 14b. In the zone of matrix cracking, the intact stuffers act roughly as 
bridging springs, so p(u) will be an increasing function. When p = p c , a value close to the ultimate 
strength for loading along the stuffers, stuffer failure will commence. The corresponding value u e 
of u will be relatively small, because stuffers deform elastically for u < u e . At higher 
displacements, p(u) will fall monotonically, vanishing when u = u c at the trailing end of the zone of 
tow pullout. 


In describing the extension of a single band of damage from a stress concentrator, the brief 
rising part of p(u) has no significant role. (Material properties that do depend on the rising part of 

p(u) are reviewed in [20]). Since u c » u e for the composites studied here, p(u) can be idealized by 
the solid line in Fig. 14b. 

The fracture process could also in principle be influenced by the work required to advance 
the zone of matrix cracking, i.e., to form new interfiller cracks. However, this energy is much 
smaller than that required to pull out failed stuffers. This suggests the closure condition that no 
stress singularity should exist at the tip of the zone of matrix cracking. The much studied cohesive 
zone model with tension softening cohesive forces then results. 

In such cohesive zone models, a fundamental role is played by the characteristic length, 
/ C h, defined by [21-23] 


The length of the cohesive zone is given to order magnitude by / c h- When the specimen width is 
less than / c h, the material will be notch insensitive, with the ultimate strength approximated by the 
condition that the net section stress should equal p c (e.g., [23,24]). In large specimens, notch 
sensitivity exists only for notches much larger than / c h. 

Estimates of E c , u c , and p c for the subject woven composites and the resulting values of / c h 
are shown in Table 10. The values of Ec are taken directly from Table 4. The estimates of p c are the 
strengths of Table 4. The ranges shown for u c correspond to the range of measured tow pullout 
lengths. The resulting values of / c h are substantial fractions of 1 m. If such high values are 
realized, then in all practical applications, these materials will be notch insensitive and damage 
tolerant under tensile loads. 


Table 10 


Estimates of Cohesive Zone Lengths for Tension Cracks 


Lightly Compacted Composites 
Heavily Compacted Composites 


Ec 

(GPa) 

30 

80 


Uc 

(mm) 

1-5 

1-5 


Pc 

(GPa) 

0.3 

0.9 


fch 

(cm) 

10-50 

10-50 


The length / c h will be reduced if the effective pullout length is reduced. The pullout length 
will decrease if the distribution of flaw strengths decreases in width, vanishing when all tows fail 
at a single strength [18,19]. If tow rupture and sliding are accompanied by disintegration of the 
surrounding composite, bridging tractions will be lost at smaller tow displacements and the 
effective pullout length will again be reduced. 

If u c is reduced by a factor of 10 from the lower end of the estimated ranges of Table 10, 
then significant notch sensitivity would arise in practical applications. 

On the other hand, damage tolerance in tension can be enhanced by appropriate composite 
design. The critical displacement u c and therefore / c h can be made large by 1) using coarse 
reinforcement (large tow denier) and 2) building in spatially distributed flaws in aligned tows. In 
the latter regard, the conjecture that tow tensile strength is lowered by lateral loads (which can be 
imposed intentionally through the choice of weave architecture) is worth pursuing. Rarely has a 
class of materials promised such scope for tailored properties. 
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Compression 


The high compressive strains to failure measured in lightly compacted composites contrast 
with the brittleness of stitched quasi-isotropic laminates [9]. The principal distinction between the 
two groups of materials lies in the distribution of flaws. Stitched laminates are heavily compacted 
and the laminas in them are relatively flat and well aligned. The critical stress for kink band 
formation is therefore relatively high and such geometrical flaws as exist are distributed narrowly 
in strength. Consequently, the first kink failure occurs usually in an outside ply at a stress near the 
critical stress for kinking of its neighbor. The stress concentration around the first failed ply causes 
immediate failure of the neighbor; and so the damage propagates unstably across the specimen. In 
the lightly compacted 3D woven composites, geometrical flaws are distributed broadly both 
spatially and in strength. The first tow failure is generally associated with a severe flaw. The 
resulting stress concentration may not suffice to fail neighboring tows, which may have no major 
flaws in that vicinity. The next failure event is then at a distant site uncorrelated with the first; and 
so damage accrues in a spatially distributed fashion and the material is ductile. 

The heavily compacted composites are appreciably more regular than the lightly compacted 
composites. Heavy compaction has minimized distortion of stuffers, the principal load bearing 
tows. However, residual distortions such as those depicted in Fig. 2 still constitute geometrical 
flaws more severe than those found in stitched laminates. Consequently, one would expect the 
heavily compacted composites to display damage tolerance and ductility in compression lying 
between those of the lightly compacted composites and stitched laminates. This is indeed the case. 

Compressive strength, on the other hand, is likely to be degraded by geometrical 
irregularity. Table 1 1 compares the ratio of compressive strength to fiber volume fraction for the 
lightly and heavily compacted composites. In composites of equal irregularity, this ratio would be 
expected from Eqs. (4) and (5) to be constant. It tends to be much greater for heavily compacted 
than for lightly compacted composites. 

Thus the lightly and heavily compacted composites demonstrate an essential balance in the 
properties of woven composites. While reduced geometrical irregularity in the latter leads to higher 
compressive strengths than would follow from consideration of volume fraction alone, this comes 
at the price of reduced strain to failure in compression. 
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Table 11 


Variation of Compressive Strength with Fiber Volume Fraction 



Strength a /f s V b 


Strength a /f s V b 

Preform Label 

(GPa) 

Preform Label 

(GPa) 

/-L-l c 

1.6 

h- L-l 

1.8 

/-L-2 C 

1.0 

h- L-2 

2.1 

/-T-l c 

0.9 

h- T-l 

1.6 

/-T-2 C 

1.1 

h-T-2 

1.6 

1-0 

- 

h-O-l 

1.8 



h-O- 2 

1.7 


a Strength for dog-bone specimens. 
b Product f s V from Table 2. 
c Strength from Table 4. 

Apart from questions of flaw distribution, compressive failure in the heavily compacted 
composites is also affected by the occasional observation of delamination and subsequent Euler 
buckling. This failure mode is not observed in either the lightly compacted composites or in 
stitched laminates, where any delamination cracks that occur remain small rather than propagating 
the length of the specimen and there is no buckling of the material on either side of them. A simple 
model of the role of through-thickness reinforcement in bridging such delamination cracks has 
been presented in [25] and [26]. Under uniaxial loading, buckling must precede delamination crack 
growth, since otherwise the specimen remains flat and there is no crack driving force. Through- 
thickness reinforcement raises the critical stress for buckling by providing an elastic foundation for 
the material on either side of any delamination crack. If buckling is thus suppressed, the 
delamination crack will not propagate and failure proceeds via kink band formation. 

The efficacy of the through-thickness reinforcement depends on the effective stiffness it 
supplies in acting as an elastic foundation, which depends on the denier and fiber type of the 
through-thickness tows and their volume fraction. A lower bound to the foundation stiffness can 
be calculated by assuming that the through-thickness tows transfer load to the rest of the composite 
only at their extrema. In stitched laminates or through-the-thickness weaves (preforms whose 
labels contain T or O in Table 1), the extrema lie at the specimen surfaces. Simple estimates can 
then be made of the minimum volume fraction required to suppress buckling at loads up to that for 
failure by kink band formation. This critical volume fraction is independent of the delamination 
length [25,26]. For current stitched laminates, the stitching density is probably one or more orders 
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of magnitude greater than required to meet this criterion [26], In this regard, the stitched laminates 
are overdesigned. 

In woven composites with the volume fractions of Table 2, the same ought to be true [26]. 
However, the warp weavers in woven composites are distorted during consolidation, especially, as 
one would expect, in the heavily compacted composites (Fig. 4). The distortion can greatly reduce 
the effective stiffness of the warp weavers under through-thickness loads: they deform in shear 
rather than in tension. The degradation is evidently so severe in the heavily compacted composites 
that the elastic foundation can become too soft to suppress buckling, and failure by kink band 
formation is sometimes not achieved. Confirmation of this conjecture was sought by measuring the 
longest delamination crack visible on the sides of specimens subjected to approximately the same 
applied compressive strain. These measurements are plotted in Fig. 15 against the ratio of the net 
composite thickness (Table 2) to the average thickness of the dry fiber preform as received from 
the weaver. The severity of warp weaver crimp seen on sections of specimens is in qualitative 
correspondence with this ratio. So too is the propensity for delamination, as evident in Fig. 15. As 
expected from differences in constraint. Fig. 15 also shows a tendency for delamination to be more 
pronounced in dog-bone specimens than in cuboidal specimens. Figure 15 clearly demonstrates the 
virtue of achieving high fiber volume fraction by controlling the weaving process rather than 
applying heavy compaction loads during consolidation. 
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Fig. 15 Maximum delamination crack length observed as a function of the degree of through- 
thickness compaction achieved in consolidation. The preform type and compressive 
strains at which measurements were made are marked at the top of the figure. 


Compression After Impact 

Provided failure occurs by kinking rather than by delamination and buckling, a simple 
extension of the idea of the cohesive zone used for notch sensitivity in tension could be used to 
estimate residual compressive strength after impact (e.g., [28]). 

Consider the impact damage zone to be a hole, which is an extreme assumption of the 
damage state and should lead to conservative strength estimates. Upon compressive loading, 
kinking will propagate out from this stress concentrator, much as tow rupture would progress in 
tension (Fig. 14a). Bridging tractions, in this case propping the crack open and shielding 
unkinked stuffers ahead of the damage zone from compressive loads, will arise because of the kink 
lock-up mechanism discussed below in Section 6 and Appendix D. This amounts to a compressive 
cohesive zone, which will have some characteristic length / c h in analogy to that of Eq. (11). We 
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have insufficient knowledge of the mechanics of kink initiation and stress redistribution following 
kinking to make good estimates of / ch . However, our qualitative observations of damage 
mechanisms suggest that compressive cohesive zones of the order of several cm. are feasible for 
these coarse woven composites. In that case, the success of net section considerations in estimating 
strength after impact in Section 3 would be no surprise. The damage zone there was only = 1 cm. 

Further research on this idea would seem warranted. 

4. Fatigue 

4.1 Compression-Compression 

Strain-life tests were executed under uniaxial, compression-compression loading using 
dog-bone specimens. Results for lightly compacted composites are shown in Fig. 16, with the 
data for one cycle being the monotonic data of Table 4. 



Cycles to failure 


Fig. 16 Stress-life data for lightly compacted woven composites under uniaxial compression- 
compression fatigue. 

Woven composites and stitched laminates have a similar contrast in fatigue to that observed 
in monotonic loading. A few fatigue tests for stitched laminates showed very flat S-N curves: an 
increase in load of 10% (from 38 ksi to 42 ksi for the specimens tested) reduced a fatigue life 
exceeding one million cycles (run-out in test) to failure in one cycle. Figure 16, in contrast, shows 
significant slope to the S-N curves. 
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Figure 17 shows stress-life data for four of the high volume fraction (heavily compacted) 
woven composites supplied by Boeing. The tests show that the orthogonal interlock architecture 
yields the highest strength, while the layer-to-layer angle interlock architecture yields the lowest 
strength. This is consistent with our view that angled warp weavers act as geometrical flaws where 
they wrap around fillers, which can be quantified by measurements of misalignment angles for 
stuffers. The layer-to-layer angle interlock architecture has the highest misalignment angles; the 
orthogonal interlock architecture the lowest. However, there is also substantial delamination in 
some of these weaves during fatigue, apparently depending on the extent to which warp weavers 
have been crimped during processing. Crimping reduces the resistance to delamination by lowering 
bridging tractions across incipient delamination cracks. 



Fig. 17 Load-life data for heavily compacted composites under compression-compression loading 
(load ratio R = -<»). 

The mechanics of failure in fatigue and monotonic compression. are evidently different. 
Under monotonic loading, stuffer failure was preceded by extensive debonding of both warp 
weavers and stuffers from the surrounding matrix, with the former debonding events apparently 
initiating the latter in many cases. In fatigue, the order of events was quite different. Figure 18 
shows the state of damage observed on the gauge section in a micrograph that was fortuitously 
shot just a few cycles before failure in a specimen of type £- L-l. While failed stuffers can easily be 
seen, all warp weavers remain well bonded to the surrounding matrix. There is indeed very little 
matrix cracking, in striking contrast to observations for monotonic loading. This feature has been 
confirmed by sectioning fatigued specimens. While sectioning revealed interior stuffer kink band 
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Fig. 1 8 One of the cut surfaces of a specimen of layer-to-layer angle interlock weave (^-L- 1 ) just 
prior to failure in compression-compression fatigue. 

failures, the only matrix cracking appeared to be caused by those stuffer failures, rather than being 
a cause of stuffer failures, as in monotonic loading. 

We infer that stuffer failure occurs in fatigue by the accumulation of damage within an 
individual stuffer. The damage is not triggered by events external to the stuffers, such as 
debonding or other matrix cracking that might allow buckling. Instead, cyclic loading slowly 
weakens the resin within the stuffer, reducing load transfer between neighboring fibers and 
therefore the critical load for kink band formation. 

One fundamental characteristic pervades observations for both monotonic and cyclic 
loading. In both cases, woven composites exhibit damage tolerance, showing high strain to failure 
in monotonic loading and significantly sloping S-N curves in fatigue. Furthermore, post-mortem 
observations in both cases reveal damage distributed throughout the gauge section, in contrast to 
the case for stitched laminates, which fail by highly localized kink bands in either monotonic 
loading or fatigue. 

4.2 Tension-Tension and Tension-Compression 

The progression of damage under tension-tension or fully reversed loading is distinguished 
from that in compression-compression fatigue by the relatively early appearance of matrix cracks 
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normal to the applied load axis. These microcracks increase in density with cycles and occur at 
tensile amplitudes in the elastic regime of monotonic tension data. 

Although the matrix cracks cause a modest increase in specimen compliance, ultimate 
failure is the result of tow rupture. Neither is there any evidence that the matrix cracking directly 
induces tow failure. Instead, it would appear that just as in compression-compression fatigue, tow 
rupture in tension signals the maturation of degradation that is internal to the load bearing tows. 
There is again correlation between sites of rupture and misalignment. 

5. Models Appropriate to Different Aspects of Composite Behavior 

The experiments reported above and elsewhere suggest that the mechanical properties of 3D 
composites can be separated into two categories: those that can be predicted to within experimental 
scatter by elementary models; and those than cannot, i.e., those for which our Binary Model is 
designed. Table 12 summarizes this classification. 

Problems for Which the Binary Model is not Required 

Flat panels consisting of a single weave type behave as orthotropic bodies in the elastic 
regime. Elastic constants are well approximated by combining rules of mixtures or other simple 
models of unidirectional composites with orientation averaging models (or even standard laminate 
theory) and some crude estimates of the softening effects of tow irregularity. 

Unnotched strength in tension can be estimated from the strength quoted for pristine fibers 
by the manufacturer, corrected for volume fractions. Unnotched strength in compression is 
governed by the mechanics of kink band formation, following Eq. (4). 

The extent of fiber pullout observed in unnotched tension tests (e.g., Fig. 10) suggests that 
notch sensitivity in tension should be modeled via a cohesive zone of damage extending from any 
stress concentrator. Within the cohesive zone, the mechanics of tow pullout will govern the relation 
p(u) between the bridging tractions, p, acting across the damage zone and the displacement 
discontinuity, 2u (or crack opening displacement, u). In Section 3, the characteristic cohesive zone 
length, / c h, was shown to be to order of magnitude 0. 1-0.5 m. The function p(u) can be deduced 
directly from measurements of force and displacement on unnotched tensile specimens, as long as 
the specimen width is much less than / c h, e.g., ~ 10 mm. Once p(u) is known, notch sensitivity 
and the influence of part size and geometry on strength can be computed from the relatively simple 
and well developed fracture mechanics of cohesive zones or bridged cracks. 
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Table 12 


Predicting the Properties of 3D Composites 

(i) Some Properties Predicted by Simple Models 


Property 

Model 

Stiffness of flat coupons 

Rule of mixtures/mean field models 
Laminate theory 

Rough estimates of the effects of tow irregularity 

Unnotched strength for aligned loads 

Compression: criterion for kink band formation 
Tension: tow rupture strength 

Notch sensitivity/fracture toughness in 
tension 

Cohesive zone model 

Delamination and buckling 

Beams or plates on an elastic foundation 


(ii) Some Problems Requiring a Computational Model (the Binary Model) 


Problem 

Remarks 

Stiffness/strength of integral structures 

Require stress distribution in tows in complicated 
arrangements 

Progression of damage in monotonic 
loading and fatigue 
Localization/delocalization of damage 
Open-hole compression 
Fatigue near stress concentrators 
Constitutive law for cohesive zone in 
tension 

Depend on local stress distributions, distributions of 
flaws, and load redistribution following local 
failure — stochastic problem 


Delamination and subsequent buckling under monotonic compression, the primary 
mechanism of failure against which 3D reinforcement has been introduced into polymer 
composites, can also be modeled relatively simply. The delamination problem can be modeled by a 
variant of existing laminate theories in which elastic springs couple separable laminae (e.g., [29]). 
For 3D composites, the springs represent the through-thickness reinforcement. The problem of 
buckling of delaminated layers in 3D composites can be described as that of classical buckling 
plates on an elastic foundation [25,26]. 

For nearly all of the above properties, the 3D composite behaves essentially as a laminate of 
homogeneous layers. The effects of the 3D weave or of irregularity in tow positioning are either 
small or are determined by averages over large volumes of material. The sole exception is modeling 
based on the concept of a cohesive zone, for which the crucial relation p(u) depends strongly on 
the irregularity and geometrical details of the reinforcement. However, engineering predictions 
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based on the cohesive zone model can be completed by determine p(u) experimentally; for this 
particular purpose, the micromechanics underlying p(u) need not be modeled in detail (e.g., [30]). 

Problems to be Solved by the Binary Model 

Other important problems defy such relatively simple modeling. They are generally those in 
which macroscopic behavior depends on the details of load distribution throughout the composite. 
Some examples are given in part (ii) of Table 12. 

One very important application of 3D composites is the fabrication of integral structures. 
Two examples from weaving technology are integral box beams, containing predominantly axial 
yams (stuffers) in the upper and lower surfaces, with ±45° yams (warp weavers) in the sides 1 ; and 
integrally woven skin/stiffener panels (e.g., [31]). Because of the complex reinforcement 
architecture in such structures, orientation averaging models are unlikely to be reliable even in the 
elastic regime, especially at critical junction regions such as where two sides of a box beam or a 
stiffener and skin merge. Predictions of stiffness and strength require the calculation of loads in 
geometrically complex arrangements of tows. 

While the effects of complex tow arrangements present a deterministic problem, several 
other aspects of composite behavior depend on how random flaws are distributed in both strength 
and space. These properties demand a stochastic model. 

In compression, the more misaligned segments of tows will fail by kink band formation at 
lower values of the local axial stress, following Eq. (4). If lateral loads also act on a tow as from 
tow wrap-around, an additional shear stress, T/, is induced. For the simplest assumptions 
concerning local fiber and kink band geometry, a k is lowered further according to [15] 

Ok = K - 1/] / <t>k ^ 

The tensile rupture of tows has not been quantitatively modeled. However, it seems 
plausible that segments of tows that are unusually bent, squashed, or subjected to lateral loads will 
have reduced tensile strength. 

In both compression and tension, fluctuations in tow alignment also affect the onset of 
nonlinearity by causing uneven load distribution. A tow segment with unusually high waviness is 


'Ray Edgeson, Cambridge Consultants, Cambridge, UK; private communication, 1992. 
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more compliant under axial loading than one that is already straight. Straighter tow segments will 
therefore bear a disproportionate share of the load and tend to fail early in loading. 

Thus the onset of damage, whether under monotonic or cyclic loading, must depend on the 
distribution of flaw strengths and the evenness of load distribution. The progression of damage at 
higher strains will depend on both of these factors as well as the way loads are redistributed around 
a local failure event. In both compression and tension, load redistribution around a kinked or 
ruptured tow is mediated by friction acting around the periphery of the broken tow. The critical 
stress for frictional sliding dictates the distance along the failed tow over which the tow is reloaded 
to far field loads by load transfer, and it therefore dictates the range of interaction of flaws. 

One important characteristic of damage progression is whether successive local failure 
events form a localized band of macroscopic damage or whether they are delocalized and widely 
distributed over the gauge section. A transition from localized to delocalized damage, manifested as 
a brittle-ductile transition in compressive stress-strain curves, has already been noted for 3D 
stitched and woven composites (Section 3). Quantitative analysis of the transition requires the 
binary model, with an appropriately detailed solution of the statistics of local failure events and 
load redistribution. Modeling the transition from localized to delocalized damage is the key to 
modeling strain to failure and damage tolerance. 

While predictions of notch sensitivity in tension can be made using an empirically 
determined constitutive law, p(u), in a cohesive zone model, composite design requires 
understanding how p(u) is determined by microstructure. Tow pullout lengths are determined 
partly by the flaw distribution within tows, wider distributions favoring long pullout lengths (e.g., 
[18,19]. Pertinent flaws comprise both intrinsic strength variations and geometrical irregularities, 
especially locations where cross tows might impose weakening lateral loads on an aligned tow. 
Pullout loads are also influenced by transverse compression experienced by ruptured tows, since 
pullout is resisted by friction. Transverse compression can be strongly enhanced by through- 
thickness tows, which often survive the rupture of neighboring axial tows, by the mechanism 
illustrated in Fig. 19. Computing how all these factors of reinforcement architecture and 
irregularity influence p(u) is another application of the Binaiy Model. 

The likely success of cohesive zone models for predicting notch sensitivity in tension is due 
to the large values estimated for the characteristic cohesive zone length, / ch . Since damage is spread 
over such large lengths before failure, details of the tow geometry near the stress concentration 
have minimal effect. Stochastic quantities are sensed only in their averages over the cohesive zone. 
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Fig. 19 Schematic of the development of transverse compression by warp weavers during tensile 
loading. White arrows indicate the sense of loads imposed on neighboring tows by warp 
weavers, which tend to straighten when strained in the direction of the load, even in the 
elastic regime. The case shown is a layer-to-layer angle interlock weave. 

In compressive loading or cyclic loading, where tow pullout lengths may not be as large as in 
monotonic tension, the failure process could depend more critically on events confined to a 
relatively small volume near a stress concentrator — perhaps containing only a few tows. If 
experiments show such behavior (they are in progress), the Binary Model will be required to 
compute the effects of notch shape, tow positioning, and random geometrical flaws. 

6. Binary Model of a 3D Composite 

The Binary Model is a finite element model, in which the highly anisotropic and 
heterogeneous structure of a 3D composite is resolved into simple constituents: reinforcing tows , 
which primarily represent the axial properties of individual tows; and an “effective medium , 
which represents all other properties of the tows, resin pockets, voids, etc. in an average sense. 
The usefulness of this division rests on the fact that the axial modulus of the reinforcing fibers, 
whether graphite, glass, or other material, is generally two orders of magnitude greater than the 
modulus of the resin. The axial modulus of the fibers dominates strength and stiffness under 
aligned loads, while the modulus of the resin dominates properties that depend on the effective 
medium, notably the shear and transverse stiffnesses and Poisson’s effect. When the model is 
discretized, the tows are divided into two-noded line elements possessing axial rigidity only, with 
no prescribed shear or bending resistance. The effective medium is divided into solid elements, 
which, at least in the elastic regime, are defined to be homogeneous and isotropic. The effective 
medium elements and the tow elements are coupled by imposing constraints between certain nodes 
of each. The constraint will usually comprise an undamaged state, in which the nodes simply share 
the same coordinates, and a damaged state, allowing some relative displacement. As in the real 
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composite, no two reinforcing tows are coupled directly. They interact only via the effective 
medium. Tow and effective medium elements are commensurate with the characteristic scale of the 
reinforcement architecture, e.g., the distance between points at which one tow crosses two other 
tows successively. Relatively large sections of the composite structure can thus be modeled in a 
calculation with a modest number of degrees of freedom. In most cases, the volume modeled is 
much greater than that of the hypothetical unit cell from which a geometrically ideal model could be 
constructed by translation and repetition. In applications to date, the size of each model calculation 
has usually been determined by the gauge section of some specimen being simulated. 

Both tow and effective medium elements are nonlinear, with plasticity and local failure 
incorporated in their assigned constitutive properties. The constitutive laws for tow and effective 
medium elements also embody stochastic parameters. 

Figure 20 shows a typical arrangement of nodes on tow and effective medium elements in 
a small volume of a layer-to-layer angle interlock woven composite. In this particular architecture, 
stuffers and fillers lie in orthogonal layers, while warp weavers supply through-thickness 
reinforcement by looping above and below individual fillers in adjacent filler layers. This structure 
lends itself to the cuboidal effective medium elements exemplified by the shaded volume in 
Fig. 20. 

As in any discretization of a continuous (or piece-wise continuous) body, there is some 
arbitrariness in the choice of element size. The choice illustrated in Fig. 20 entails the minimum 
density of tow nodes required to reproduce the topology of the reinforcement faithfully. A higher 
density of nodes could be chosen, but that would betray the spirit of finding the simplest possible 
realistic formulation. Nevertheless, as well as the assurance of computational precision in modeling 
tows, other physical conditions must also be met in choosing the element size. These will be 
discussed below. 

Once the element size has been chosen, the length scale it introduces determines the way 
geometrical fluctuations are treated. Thus, the irregularity of a tow that undulates over wavelengths 
larger than the element size can be mimicked simply by displacing appropriate nodes on that tow in 
the initial, load-free configuration of the model. If a geometrical property fluctuates over a length 
less than the element size, the variation must be incorporated in the constitutive law provided for 
that element. Some explicit examples of this will be given below. 
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Fig. 20 Tow and matrix elements in a layer-to-layer angle interlock woven composite. 

The conceptual division of the composite into tow and effective medium elements is 
ultimately a matter of convenience. The end product is a nonlinear finite element model, which 
must be equivalent to a model containing only solid elements, with the effects of the tow elements 
incorporated in appropriate anisotropic and nonlinear constitutive laws. However, the geometrical 
complexity of the reinforcement would oblige the definition of many different types of solid 
element in such a model, depending on the local tow configuration. Furthermore, the treatment of 
stress redistribution around a failed tow and the computation of local axial stresses in tows would 
be cumbersome. 

As in any finite element model, the boundary conditions remain unspecified in the general 
formulation. In each application, either load or displacement conditions are specified for boundary 
nodes. There is no restriction in principle on the distribution of the boundary forces or 
displacements; or on the shape of the boundary itself. 
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6.1 Constitutive Laws for Tow Elements 

The properties of tow elements follow from various elementary arguments. 

Elastic Properties 

The axial elastic modulus, E t , of a tow that is initially straight can be estimated by the rule 
of mixtures: 


E t = V t E f + (l-V t )E r , (13) 

where V t is the volume fraction of fibers within a single tow and Ef and Ef are Young’s moduli for 
the fibers and resin, the former measured axially. The axial stiffness, k t , of the corresponding tow 
element in the binary model is given by 

= (Et - E m ) At (14) 

where A t is the tow’s cross-sectional area and E m is Young’s modulus for the effective medium 
(specified below). The subtraction of E m in defining kt avoids double counting that would arise 
because the effective medium elements fill all space. The area A, is deducible from the length per 
unit mass, y, of the tow (known as the “yield); the density, pf, of the fibers; and Vt according to 

Al = v^ • 05) 


If the tow is undulating initially, its response to axial loads a t can be described by the 
differential axial displacement, u, between two points separated by some gauge length L: 

u = Lc t /E t + Le u = Le 

where £„ is a strain contribution arising from straightening of the undulations. If the tow undulates 
over periods greater than the element size, both of these terms will be computed as part of the 
solution of the discretized model, with tow undulations entered explicitly as initial nodal offsets. 
However, if the tow undulates over periods less than the element size, the second term in Eq. (16) 
must be computed in advance by micromechanical modeling. The tow elements then possess a 
possibly strain-dependent, reduced effective axial stiffness, E e = da t /de. 

Calculating E e accurately is difficult. However, knockdowns due to waviness or 
irregularity should never lower the composite modulus by more than 10-20%: if they do, the 
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composite is either poorly designed or badly fabricated. Experimentally measured scatter in the 
macroscopic composite modulus is usually 5-10% in 3D woven or braided composites, while tow 
waviness itself is difficult to characterize and its statistics difficult to measure. Therefore, it is 
unlikely to be profitable to employ more than rough estimates for Ee- More usefully, simple models 
of the effects of tow waviness or irregularity allow competing composite designs to be compared 
or the effects of irregularity induced by processing to be estimated. Thus, E e is simply written 


where x is a random variable for the tow elements in a single simulation and typically x > 0.9. 

For large strains, Eq. (17) will generally be nonlinear: an undulating tow will stiffen as it 
straightens. However, in most applications, it will be valid to assume that the knockdown factor % 
is independent of strain. At strains high enough for tow straightening to change E e , the composite 
stress-strain relation is likely to be dominated by tow and matrix failures. 

In representing the resistance of the tows to lateral deflections by the stiffness of the 
effective medium, it is assumed that such deflections would arise from shear alone. However, at 
least in principle, both shear and pure bending could contribute significantly to lateral deflections. 
The proportions of the total deflection arising from each will depend on the length of the tow 
element, among other things, with long elements favoring the dominance of the pure bending 
contribution. A simple estimate of these effects is presented in Appendix C. For the computational 
tow elements defined here and for all current applications in polymer composites, shear is the 
dominant mode of lateral deflection. Thus shear stiffness is properly ascribed to the effective 
medium, leaving the tow elements themselves with no inherent resistance to lateral loads. This 
division has the added virtue of minimizing the degrees of freedom in the model. 

Strength in Compression 

The strength of a tow in compression is dictated by the mechanics of kink band formation, 
with the critical value, or, of the axial stress being given by Eq. (4). Thus a k depends on two 
factors, the misalignment angle <> k and the local shear stress due to lateral loads, X/, both of which 
may change as the composite responds to load. The shear stress, X/, in any tow element can be 
computed at each load increment from the shear stresses in adjoining effective medium elements. 

How the misalignment angle <t> k is treated depends on whether the misalignment occurs 
over a gauge length that is larger or smaller than the tow element size. If it is larger, then <> k is 
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computed from the nodal displacements. If smaller, then the role of is subsumed in the 
constitutive properties of the tow, which are prescribed a priori from micromechanical arguments. 
Misalignments can be measured from digitized photographs of specimen cross-sections: the 
inferred <f»k will generally be a random variable. 

Strength in Tension 

The stress a t in any tow element can be strongly affected by tow undulations. For 
undulations of wavelength greater than the element length, resulting variations in at along the 
length of a tow are computed directly in solving the model. Undulations lying within a single 
element are represented by the reduced effective stiffness, Ec, of Eq. (17): thus 

<*t = E e et, ( 18 ) 

where £t is the axial strain in the tow element implied by its nodal displacements alone. 

The failure of a tow in tension is modeled simply by the criterion that failure occurs when 

a t“ CT t • (19) 

where at is the axial stress in the tow element and the critical stress, a[ c \ is a material property. 
The critical stress depends upon intrinsic flaws in the tow, including flaws associated with 
crimping or distortion caused by consolidation. It is a random variable whose distribution of initial 
values will probably always be evaluated by fitting the model to experimental data. It may be 
reduced during loading by lateral loads imposed on a nominally aligned tow by neighboring tows. 
If so, the value of a| c) for elements in the aligned tow could be lowered during a simulation in 

proportion to such lateral loads, as measured by the maximum shear in adjacent effective medium 
elements. 

Post-failure Properties 

The most important local phenomenon following failure of a tow in either compression or 
tension is the transfer of load to neighboring tows. This is represented by the constitutive law 
coupling tow and effective medium nodes, as described in the following section. Once the tow has 
failed, the stiffness of the failed element is usually reduced to zero. 
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However, there are circumstances where experimental evidence indicates more complex 
behavior. One is intimated by the complex kink band structure of Fig. 21. Kink bands begin 
forming at strains of 1-2%, whereas ultimate failure in tests occurs at strains up to 15%, 
whereupon post mortem examination is undertaken. Multiple kinks, such as those in Fig. 21, have 
presumably developed as a succession of kinking events over such large strains. One explanation 
of this phenomenon is that kinks lock up after a certain amount of axial strain, whereupon the tow 
can again bear large loads. Further loading can lead to a new kink band, which is very likely to 
abut the prior damage, since that will act as a nucleation site. The mechanics of lockup and its 
relation to axial sliding of a failed tow are described in Appendix D. The upshot is the schematic 
stress-strain response for the tow element shown in Fig. 22. 
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Fig. 21 System of kink bands at one site of local failure in a woven interlock composite tested to 
failure in uniaxial compression. 
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Fig. 22 Schematic of the stress-strain history of a tow element within which multiple kink 
banding occurs in compression. 

Properties in Fatigue 

The observed absence of microcracking prior to kink band formation in compression- 
compression fatigue (Section 4) implies that fatigue damage in compression consists of 
degradation of the interior of the primary load bearing tows. It has previously been pointed out 
[33] that fiber misalignment causes lateral or resolved shear loads in unidirectional composites that 
are large enough to damage the epoxy resin. In compressive fatigue loading, such damage will 
gradually lower the effective shear strength of the resin, which in turn will lower the critical stress 
for kink band failure (Eq. (4)). When that stress falls below the maximum applied compressive 
stress, the tow will fail. This model is consistent with the observation that kink bands form earliest 
in fatigue at locations of maximum misalignment. A feasible fatigue law for compressive fatigue is 
that 


^J- = -A,(Aa 1 <t) k ) n ' (A, > 0) , (20) 

where N is the number of elapsed fatigue cycles, Aa t the local stress amplitude in a tow element 
whose misalignment is <f>k, and Aj and ni are to be evaluated empirically. The exponent nj could 
be deduced from the slope of a strain-life curve by varying the degree of misalignment achieved in 
processing or by comparing the cycles to kink band formation at sites that differ in misalignment. 
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In tensile fatigue loading, the first evident fatigue damage is matrix cracking normal to the 
load axis (Section 2.2). This will be dealt with below under the constitutive properties of the 
effective medium. Matrix damage is followed by tow rupture. The rule of mixtures demonstrates 
that the increase in load in the aligned tows because of matrix cracking must be slight. More 
importantly, tow misalignment might be expected to introduce significant shear stresses within 
tows in tension just as in compression, leading to direct fatigue damage of the resin within tows 
[33]. This might lower tow strength if fibers suffer attrition following fragmentation of the resin. A 
law similar to that in compression is suggested: 

J (c) 

^=-A 2 (AG t <l> k ) n2 (A 2 >0), (21) 

dN 

where A 2 and n 2 are further empirical parameters. 

Whether Eqs. (20) and (21) are confirmed by experiments or whether tests will suggest 
alternative forms remains a topic of research. 

6.2 Constitutive Laws for Effective Medium Elements 

Elastic Properties 

An assembly of tow elements such as that of Fig. 20 cannot adequately model macroscopic 
shear stiffness, through-thickness stiffness, or Poisson’s effect. The effective medium corrects 
these deficiencies. 

One face of a typical effective medium element is outlined on the micrograph of an angle 
interlock composite shown in Fig. 23. It contains resin pockets and parts of tows oriented in 
various directions. While its elastic properties are complex in detail, those that remain after the axial 
stiffness of tows has been removed to tow elements can be approximated very simply in their 
spatial average. For most composite properties, it is a fair approximation to assume that the 
effective medium is homogeneous and isotropic in the elastic regime, with properties given by 
rules of mixtures. Let G f and G r be the shear moduli and v f and v r Poisson’s ratios for the fibers 
and resin. Then the shear modulus G m and Poisson’s ratio v m for the effective medium can be 

written 

_L_ Vf LVf (22a) 

G m -G f G r 
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and 


v m ~VfV f + (l - V f )v r , (22b) 

where, making due allowance for resin pockets and fluctuations in tow density, Vf is the volume 
fraction of all fibers averaged over the composite, as measured, for example, by weighing the 
fibers after removing the resin by acid digestion. 
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Fig. 23 A section of a through-the thickness angle interlock woven composite, consisting of AS4 
carbon tows in an epoxy resin matrix. The rectangle outlines one face of a typical 
effective medium element. ~ J 

Strength 

Tensile failure of effective medium elements occurs, for example, in the delamination of 
layers of stuffers and fillers during compression or in the formation of matrix cracks normal to the 
load axis under monotonic or cyclic tensile loading. The strength of an effective medium element is 
denoted O^. Examination of many micrographs reveals that the most important cracks form in 
layers of resin between tows rather than within the tows themselves. 2 Therefore, is likely to 
reflect the properties of the resin. However, will also depend strongly on the geometrical 


2 Some microcracks are formed within tows when high through-thickness compaction loads are applied during 
processing. Such cracks have no effect on damage progression that has been observed as yet, although they could 
conceivably lower the critical stress for kink band formation and accelerate fatigue damage within tows. 
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details of local fluctuations in tow deployment, resin porosity, etc., which are extremely difficult to 
measure or model. Furthermore, the grid is so coarse that the value of at which matrix crack 
propagation occurs for a given value of the applied load is likely to depend on the element size. 
Therefore, will usually be treated as a model-dependent, empirical parameter. 

In fatigue, one might conjecture a law paralleling Eqs. (20) and (21): 

^ = -A 3 (Ac m ) n3 (A 3 >0), (23) 

dN 

where A 3 and n 3 are empirical parameters and A<r m is the cyclic stress amplitude in the effective 
medium element. 

For failure in either monotonic loading or fatigue, experiments show that the fracture plane 
almost always either separates pairs of adjacent tows or separates layers of tows. For a model 
geometry such as that of Fig. 20, the relevant component of stress to be compared with (jj„ will 
accordingly lie in one of the Cartesian directions shown in Fig. 20. It should be averaged over the 
tow element, since computed variations within an element depend on the choice of element size. 

Post-failure Properties 

After failure, an effective medium element will have no remanent strength in tension, but 
will continue to support load in compression. It can also bear tensile loads in directions orthogonal 
to the plane in which it failed. For example, the microcracks observed normal to the load axis in 
tension-tension fatigue diminish the axial stiffness but do not necessarily imply the delamination of 
stuffers and fillers. Thus, after failure, effective medium elements are anisotropic. 

6.3 Constitutive Laws for Coupling Springs 

Coupling Between Tow and Effective Medium Elements 

When a tow fails, whether in axial compression or tension, stress redistribution is 
governed by sliding of the broken tow parallel to its axis in the vicinity of the failure site. 
Experimental observations suggest that sliding is Mode II displacement of a circumferential debond 
crack. A reasonable description of the redistribution of load is given by the shear lag model of 
Fig. 24(a). In the shear lag model, load is transferred from the tow to the surrounding composite 
(or “effective medium” in the binary model) via a constant frictional shear stress, Tf, acting over the 
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Fig. 24 (a) Shear lag depiction of stress redistribution near a site of tow failure, (b) Analog of (a) 
in the binary model. Tow and effective medium nodes have been drawn in (b) with 
different vertical coordinates solely to make them separately visible. Under axial loads, 
only their horizontal coordinates could differ in the orientation shown, (c) The* 
constitutive law for axial displacements coupling springs between tow and effective 
medium elements. 
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sliding boundary. The shear tractions restore the axial load, <J X , in the tow from zero at the site of 
failure to the far field 3 value, a t , over a characteristic length, / s , given by force equilibrium: 


l = 

s |dc x /dx| 

dq x | _ sif 
dx I A t 

where s is the circumference of the tow and A t its cross-sectional area. 


(24a) 

(24b) 


If the binary model is to describe stress redistribution correctly, then the gradient of axial 
stress in a tow near a failure site must have the value given by Eq. (24b), at least in its spatial 
average. This can be assured by appropriate specification of the constitutive laws for the nonlinear 
springs coupling tow and effective medium nodes. The discontinuity in the axial load in successive 
tow elements near the site of tow failure is simply the axial force, fa, imposed by the coupling 
spring. The axial load in tow elements near the site of failure is therefore the staircase function 
shown schematically in Fig. 24(b). This function will have the same average gradient as o x (x) in 

Fig. 24(a) provided 

|fal = st 0 L , < 25 ) 

where L is the computational element length. Equation (25) prescribes a force that is independent 
of the relative displacement, d a , of the relevant tow and effective medium nodes in the axial 
direction, as in Fig. 24(c). 

The relative axial displacement of the nodes it couples is the only degree of freedom needed 
for a coupling spring between a tow and the effective medium. The tow and effective medium 
nodes always coincide in their lateral displacements. 

In the common case that x 0 represents frictional sliding, its value should change with the 
transverse compression acting on the tow. The latter can be evaluated and continually updated by 
averaging the stress fields in adjacent effective medium elements during a simulation. In this way, 
the additional compressive loads introduced by warp weavers via the mechanism of Fig. 19 can be 
modeled. 

3“Far field” refers here simply to the composite beyond the domain of sliding. Since the stiffness of tow elements is 
generally a random variable, the stresses in tow elements even in domains far removed from any stress concentrator 
do not share a unique value. 
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The use of coupling springs not only provides approximately correct stress recovery in a 
tow near a site of failure, but also avoids physically improper behavior which would otherwise 
arise from the fact that tow elements have no cross section. 

Coupling Between Warp Weavers and Fillers 

In practice, the optimum combination of in-plane and through-thickness properties is 
usually attained in a 3D woven composite if the warp weavers are of considerably lighter denier 
than the stuffers or fillers and accordingly of lower volume fraction [26]. This invites the 
simplification of coupling warp weavers by springs directly to the fillers around which they wrap, 
rather than via effective medium elements, thus reducing the degrees of freedom of the model. The 
locations of such springs are shown in Fig. 20, while constitutive laws for them are given in 
Appendix E. This simplified treatment of warp weavers is satisfactory at least for in-plane 
properties in the elastic regime. In modeling failure, the critical role of warp weavers is to impose 
lateral loads on fillers, which remains well represented. 

7. Applications of the Binary Model 

The binary scheme outlined above is potentially applicable to very diverse 2D and 3D 
woven and braided composites, given suitable definitions of grids and minor modifications of the 
constitutive laws. It is also ultimately a viable approach to modeling continuously reinforced 
structures, such as integrally woven or braided skin/stiffener components for airframes. Such 

longer term goals noted, the following examples and remarks address the subject materials of this 
program. 

Monotonic Loading 

The computer code written to solve the binary model was based on the ABAQUS finite 
element packaged Some illustrative simulations carried out with the ABAQUS-based code are 
presented here. They demonstrate some of the effects of introducing random strengths for tow 
elements and randomness in tow positioning. 

The simulations were of uniaxial tension of through-the-thickness orthogonal interlock 
woven composites under displacement control. Cartoons of the model structure are shown in 
Fig. 25. Warp tow elements are shown there as ribbons, faces of effective medium elements as 
quadrilaterals, and sections of fillers as black dots. Each simulation modeled a section of material 


t Hibbitt, Karlsson, and Sorensen, Inc., Pawtucket, Rhode Island. 
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containing ten distinct layers normal to the filler direction, which Fig. 25 shows in an exploded 
view. Six of the layers contain stuffers and two contain warp weavers. The ninth and tenth layers 
contain no warp tows; they are included to avoid tow elements (which should lie along tow axes) 
being present on specimen surfaces. Simulations were executed for both ideal and irregular tow 
positioning. Irregularity was introduced by offsetting the initial, stress-free coordinates of nodes. 
The offsets were chosen by Monte Carlo methods, i.e., using a pseudo-random number generator, 
according to an ad hoc random walk model. The magnitudes of the tow displacements were on 
average about 30% of the tow spacing in any direction. 


(a) Ideal Geometry, Uniform Tow Strengths 
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(b) Random Geometry, Random Tow Strengths 



Fig. 25 Element geometry and damage sequence for orthogonal interlock composites with 
(a) ideal geometry and uniform tow element strengths; (b) random geometry and random 
tow element strengths. The inset stress strain curves contain a color code that allows the 
sequence of tow failures to be mapped. The ten parts in each exploded view are planes 
normal to the y-axis. They lie in the order indicated by the labels yi,...,yio- 
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The stiffnesses of tow elements and effective medium elements, the average tow element 
strength, and the element dimensions that were chosen are typical of AS4 carbon tows in an epoxy 
resin (Table 13). The tow elements were either assigned uniform strengths (the average (o^) of 
Table 13) or normally distributed strengths with standard deviation 20% of the average. If 
strengths were uniform, one centrally located element was assigned a slightly lower value than the 
rest, to ensure that failure started away from the specimen ends. For random strengths, the lowest 
strength value generated was always assigned to that same element. The effective medium was 
assigned infinite strength: since it is so soft, the qualitative results considered here are not greatly 
affected. 

Table 13 

Specifications for the Simulations of Fig, 23 


k s (Nt) 

2x 10 5 

kf(Nt) 

1 x 10 5 

k W w (Nt) 

2.5 x 10 4 

a x (mm)* 

2 

a v (mm)* 

1.8 

&z (mm)* 

0.65 

Em (GPa) 

10 

(a! c) >(GPa) 

2 

(°S?) (GPa) 

OO 


* See Fig. 20 for definition 


In these illustrative simulations, the degree of freedom that allows relative sliding of failed 
tow elements and the surrounding composite, as described in Section 6, was suppressed. (It is 
difficult to treat in the ABAQUS-based code.) This corresponds to the friction stress x and 
therefore the stress concentration on tows neighboring a failure site both being large (but finite). 

In each simulation, the applied strain was incremented in steps small enough that at most 
two or three and usually zero or one tow elements would fail. (Complete control is not possible in 
the ABAQUS-based code). The simulations were carried on to large strains. 

Figure 26 shows load-strain records for uniaxial tension simulations under displacement 
control in the stuffer direction. The curves show the effect of successive tow failure events. 
Figure 26(a) is the case of ideal geometry and uniform tow element strengths. Brittle behavior is 
found: when one tow element fails, propagating stress concentration causes many tow elements to 
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fail in an unstable manner, t Only the stabilizing influence of the fixed grip loading conditions and 
the unrealistically infinite strength of the effective medium prevent total failure of the specimen. 
The other three cases (Figs. 26(b)— (d)) show the effects of irregular geometry and random strength 
assignments acting separately or together. Randomness in either geometry or strength enhances 
ductility. 
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Fig. 26 


Load-strain curves for various combinations of randomness of tow strength and 
geometry. 


t In the ABAQUS-based code, all critically loaded elements must fail in a single load increment. For elements of 
uniform strength , this leads to multiple, simultaneous element failures along a single stuffer, as seen in Fig. 25(a), 
even when the load increment is very small. This unappealing and unphysical effect disappears as the distribution of 
tow element strengths broadens (Fig. 25(b)). 
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This trend is underscored by the failure sequence of tow elements, which has been 
incorporated in Fig. 25 by color coding. The value of applied strain at which any tow failed is 
revealed by matching its color against the color strips in the inset stress-strain record. 

For ideal geometry and uniform strength (Fig. 25(a)), tow failure propagates in a 
symmetric, deterministic way from the first failure site. Nearly all elements on the plane on which 
the first failure occurred fail with very little increase in applied strain. Two elements survive, 
shielded by nearby warp weavers. A complete fracture path is still formed during the first load 
drop by the failure of the elements on either side of the two survivors. The propagating effects of 
stress concentration dominate the damage evolution. 

In contrast, for the case of irregular geometry and random tow element strengths 
(Fig. 25(b)), half the tow elements on the plane of first failure remain intact throughout the entire 
process. Damage is distributed over the whole specimen, tending to occur in bursts in one distinct 
region after another. Randomness has induced a brittle/ductile transition. 

The damage sequence when only one of geometry or strength is random appears 
intermediate between Figs. 25(a) and (b), as Fig. 26 and intuition would suggest. The patterns of 
failure found for random geometry or random strength are quite similar to one another. By 
increasing the variance of either initial node positions or tow element strengths, a brittle/ductile 
transition can be induced. Which one is the stronger factor in current 3D composites remains to be 
investigated. 

Fatigue Loading 

Simulations of fatigue are similar in complexity to simulations of monotonic loading, but 
follow different constraints. Typically, either the cyclic applied load amplitude or the cyclic applied 
strain amplitude is specified and held constant throughout the simulation. Damage is then measured 
over elapsed cycles, N, treated as a continuous variable. The strength of each element decays with 
N according to Eqs. (19), (20), and (22). The parameter 4>k is a random variable assigned an initial 
value for each element. The local cyclic load amplitude, A a t or Aa m , is computed by the model. 
Equations (19), (20), and (22) are integrated to identify the first element whose strength falls in 
fatigue to a value equal to the maximum load it bears in any load cycle. That element is then failed. 
The simulation is then relaxed to determined new cyclic loads on all remaining elements, which are 
then tested for failure without further increase in N. If none fails, Eqs. (19), (20), and (22) are 
integrated again, now using the new values of local stress amplitude, until another failure is found. 
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Ultimate failure occurs when there are sufficient failed elements to cause catastrophic damage 

propagation in one cycle. Quantitative examples of fatigue simulations are being conducted in our 
new contract. 

Engineering Applications 

As illustrated above, single simulations are executed by assigning random parameters 
(strengths and irregularity) with a pseudo-random number generator. The statistics of composite 
properties are then determined from an ensemble of simulations. Ultimately, calculations would be 
more efficient in a probabilistic formulation such as a Markov chain or diffusion equation. 
However, in this investigative phase, Monte Carlo methods have the great advantage that newly 
discovered phenomena are relatively easily programmed into the solution. In particular, the number 
of independent variables can be readily changed in a Monte Carlo simulation, whereas changing the 
dimension of a probabilistic formulation is a major programming exercise. 

Input for the model consists of geometrical and material parameters, which either refer to 
deterministic quantities, e.g., average stiffnesses and lengths, or are parameters in the distributions 
of random variables. Some of the parameters can be regarded as known a priori from 
micromechanical arguments or measurements. Others will be empirical, having to be determined by 
calibrating the model against test data for the composite. A great part of current research is directed 
to determining which parameters fall into which category. 
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Appendix A 

Processing of Lightly Compacted Composites 


We consolidated the lightly compacted composites ourselves, using methods developed 
ad hoc. The matrix in all cases was formed from Tactix 138 resin and H41 hardener.^ The 
preform was placed in a reusable aluminum mould, heated to 65°C, and degassed in a vacuum of 
~ 1 Torr. The resin was then mixed with the hardener and also heated to 65°C and degassed. The 
mixture was poured over the preform (still at 65°C) and the whole assembly was degassed twice 
again. The first degassing typically resulted in bubbling out of some volatiles followed by apparent 
boiling as the pressure fell. The pressure was then cycled between 1 Torr and atmospheric pressure 
to remove small bubbles clinging to the mat. The mould was closed and the specimen cured, the 
curing cycle (chosen to maximize resin toughness) comprising 2 h at 120°C and 2 h at 177°C. All 
fabrication runs produced 25 x 10 cm panels, from which specimens were machined. 

Figure A. 1 shows representative cross sections of a through-the-thickness angle interlock 
specimen. There is no visible porosity in the optical micrograph of Fig. A. 1(a), while the scanning 
electron micrograph of Fig. A. 1(b) shows complete wetting of each fiber in an individual filler tow 
and the absence of any matrix or interfacial microcracking. If care was taken to avoid excess resin 
on the specimen surfaces (which was not especially easy, since the preform thickness varied), the 
only microcracks to be found occurred in surface pockets of resin between tows. Such microcracks 
would extend into the composite until they encountered the internal microstructure. Thus their 
average size was about the diameter of one fiber tow. The large openings of these cracks attest to 
the substantial, tensile residual stresses in the epoxy, which evidently drive the microcracking. 
However, detailed observations of failure mechanisms revealed that these initial surface 
microcracks have no role in compressive failure. 

During impregnation and cure, a small pressure was applied through the thickness of the 
preform by bolting down the lid of the mould. This pressure was sufficient to ensure firm contact 
between the mould and the thicker parts of the preform, but not sufficient to maximize fiber 
packing density elsewhere. 


§ Dow Chemical, Freeport, Texas. 
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Fig. A. 1 (a) Optical micrograph of section parallel to warp fibers through composite of preform 

/-T-l of Table 1 with Tactix 138/H41 matrix, (b) SEM micrograph of section through 
an individual filler in (a). 
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Appendix B 

Shear Flow Stress for Kinking of Tows 

This appendix addresses the shear flow stress, Tf, upon which the critical compressive load 
for kink band formation depends (Eq. (4)). 

In fact, neither the Shell 1895 nor the Tactix 138 resin alone exhibits significant plastic 
flow except under uniaxial compression. Both appear very brittle in either tension tests or “shear” 
tests using Iosipescu specimens. In the latter tests, catastrophic cracking spreads from small 
volumes of material near the notches where tensile stresses exist. Thus properties relevant to kink 
band formation in composites are not observable in the neat resin. The relevant properties are those 
of the resin when it is under the mechanical constraint of the fibers. 

One approach to measuring Tf might be to attempt shear loading experiments on tows in the 
subject 3D woven composites. However, such an experiment is hard to devise and would be even 
harder to analyze, since the local stress state is very complicated. Tests were run instead on ±45° 
laminates, following work by earlier authors [27]. 

Laminates were available for the AS4/Shell 1895 system only. They consisted of 
symmetric ±45° lay-ups of “uniweave” fabric, in which AS4 and glass fibers are combined in the 
proportion 19:1 by volume in a plane weave. “Uniweave” is prepared to facilitate handling. The 
glass fibers have an insignificant effect on failure in the present context. The total volume fraction 
of AS4 fibers is similar in the laminates to that measured within individual tows in the 3D woven 
composites. 

The laminates were loaded in uniaxial tension along the 0° direction. The onset of near- 
perfect “plasticity” in Fig. B. 1 is associated with microcracks in the resin that are arrayed linearly 
in the fiber direction (Fig. B.2). In fact, consistent with tests on neat resin, the resin near the 
microcracks shows no sign of plasticity in the continuum sense. The nonlinearity in macroscopic 
stiffness evidently arises from the microcracks themselves or frictional sliding when an array of 
microcracks collapses into a band of more severe damage at high strains. 

Since kink bands involve strains exceeding 1%, the pertinent domain in Fig. B.2 is that for 
which near-perfect plasticity prevails. The shear stress resolved in the fiber direction in any ply in 
the ±45° laminate is half the applied load. Thus, from Fig. B.l, Tf= 75 MPa for the AS4/Shell 

1895 composites. 
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Fig. B. 1 Applied load vs strain for ±45° laminates of AS4/Shell 1895 loaded in uniaxial tension 
along the 0° direction. 



Fig. B.2 A view of the original panel surface from a tested specimen, showing a linear array of 
microcracks in the resin of a ±45° AS4/1 895 laminate, the source of the “plasticity” in 
Fig. B- 1 . The array follows the local fiber orientation. 

Similar laminates of AS4/Tactix 138 were not available for testing. However, since the 
Tactix 138 and Shell 1895 resins exhibit similar properties in uniaxial tension and compression, Tf 

for Tactix 138 was also assigned the value 75 MPa. This must be a fair value for the rough 
estimates of kink band stresses in Table 7. 
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Appendix C 

Effective Flexural Rigidity of a Tow 


Consider the problem of a tow segment of length L loaded transversely by a force of 
magnitude F at one end (Fig. C.l). Suppose the other end is built in. The deflection Vb(x) arising 
from bending satisfies 

E t I^k = Fx, (C.1) 

dx 2 

where E t is the axial Young’s modulus of the tow and I is the relevant moment of inertia; while 
that, v s (x), arising from shear satisfies 

AG,^=F, (C.2) 

dx 


where G t is the shear modulus of the tow and A is its cross-sectional area. Hence the displacement 
contributions at the load point x = 0 are in the proportion 


Vs(0) _ FL / FL 3 
v b (0) AG, 3E,I 




2 


El 

G, 


(C.3) 

(C.4) 


for a tow of elliptical cross-section with semi-axis a in the direction of bending. 



Fig. C.l Cantilever beam paradigm for estimating the proportions of bending and shear in a 
transversely loaded tow segment. 
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The axial modulus E t can be approximated by Eq. (13). The shear modulus G t can be 
estimated by rule of mixtures from the shear moduli Gf and G r of the fiber and resin: 

_L_ V t 1-V, 

g, "g7 "gT ' < C5 > 

where Vt is the volume fraction of fibers in a tow. For typical graphite/epoxy or glass/epoxy 
systems, the term in Gf is negligible. 

The ratio E t /G t is typically fairly large - for AS4 graphite fibers in Shell 1895 resin, the 
material combination for many of the composites studied here, V t ~ 0.7, Ef = 250 GPa, and 
E r — 3.5 GPa, leading to Et ~ 200 GPa, Gt ~ 4 GPa, and Et/Gt ~ 50. Substituting this number 
into Eq. (C.4) shows that shear will dominate lateral deflections of a tow segment if a force couple 
acts over lengths less than five tow widths. 

This criterion is satisified for forces acting on the ends of tow elements in the binary model. 
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Appendix D 
Kink Band Lock Up 

The relation between kink band lock up and the axial sliding of debonded tow segments can 
be elucidated by simple shear lag analysis. 

Following formation of a kink band, the fiber segments within the band rotate (Fig. D.l). 
At first, the rotation is accompanied by transverse dilatation. But as rotation progresses, the fibers 
are drawn back together. At some critical rotation angle, 0 C , the volumetric strain in the band 
vanishes and further rotation is very strongly resisted. The band is effectively locked up. If the 
boundaries of the kink band form at angle p to the boundaries of the tow (Fig. D.l), the condition 
of vanishing volumetric strain leads to [31] 

0 c =2p . (D.l) 

Typically, P = 20-30°. 

Assume that kinking is accompanied by debonding of the tow over a sliding length l s and 
that the axial compressive stress at the kink band is zero prior to lock up. Then, according to the 
same shear lag model that underlies Eq. (24), sliding will produce an axial displacement u of the 
ends of each of the intact parts of the tow given by 

u = 1 °t A (D.2) 

2 E t s T 0 

where o t is the stress in the tow remote from the kink band, A t and s are the cross sectional area 
and circumference of the tow, and t 0 is the critical shear stress for sliding of the tow. 

If the length of fibers within the kink band is h (Fig. D.l), lock up will occur when 

u = k(l-cos0 c ) . (D.3) 


The corresponding value, a/, of the remote stress in the tow is 

11/2 


a/ = 


2u • E t ^ 
A t 


(D.4) 


75 



The critical stress for lock up depends via Eq. (D.3) on the kink band length, h, which is difficult 
to predict. 



Fig. D.l Schematic of fiber rotation within a kink band. 
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Appendix E 

Coupling Springs Between Warp Weavers and Fillers 

Figure E. 1 shows a schematic of a warp weaver wrapping around a filler. The warp 
weaver has a radius r w and the filler a radius rf, which can be estimated from the fiber volume 
fraction, V t , within a tow (V t ~ 70%), the denier of the tow, and the density of the fibers. In the 
elastic regime, displacement of the point P on the axis of the warp weaver relative to the point R on 
the axis of the filler is resisted by the transverse stiffness, E tr , of the two tows. This might be 
approximated by the rule of mixtures: 

E tr =[V t /Ef+(l-V t )/E r ]-«, (E- 1 ) 

where Ef and E r are the fiber and resin moduli. Assuming the contact area 4rpr w (Fig. E.l), the 
effective spring constant k w f coupling the warp weaver and filler is defined by 

k w f = E tr • 4rfr w , ( E - 2 ) 

where the spring constant relates force to proportional change in displacement. Equation (E.2) can 
be readily generalized to the case of tows containing different kinds of fibers. 

Failure of the spring occurs at some critical tensile displacement, which will usually be 
treated as an empirical parameter. 

SC-291 9-T 
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Fig. E.l Schematic of a warp weaver wrapping around a filler in an interlock weave. 
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